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This paper presents the results of a micro- 
electromechanical systems (MEMS) 
Young's modulus and step height round 
robin experiment, completed in April 
2009, which compares Young's modulus 
and step height measurement results at a 
number of laboratories. The purpose of the 
round robin was to provide data for the 
precision and bias statements of two \ 
related Semiconductor Equipment and 
Materials International (SEMI) standard 
test methods for MEMS. The technical 
basis for the test methods on Young's 
modulus and step height measurements are 
also provided in this paper. 

Using the same test method, the goal 
of the round robin was to assess the 
repeatability of measurements at one 
laboratory, by the same operator, 
with the same equipment, in the shortest 
practical period of time as well as the 
reproducibility of measurements with 
independent data sets from unique 
combinations of measurement setups and 
researchers. Both the repeatability and 
reproducibility measurements were done 
on random test structures made of the 
same homogeneous material. 

The average repeatability Young's 
modulus value (as obtained from 
resonating oxide cantilevers) was 
64.2 GPa with 95 % limits of + 10.3 % 
and an average combined standard 



uncertainty value of 3 . 1 GPa. The average 
reproducibility Young's modulus value 
was 62.8 GPa with 95 % limits of 
± 11.0% and an average combined 
standard uncertainty value of 3.0 GPa. 

The average repeatability step height 
value (for a metal2-over-polyl step from 
active area to field oxide) was 0.477 urn 
with 95 % limits of 7.9 % and an average 
combined standard uncertainty value of 
0.014 um. The average reproducibility step 
height value was 0.481 um with 95 % 
limits of ± 6.2 % and an average combined 
standard uncertainty value of 0.014 um. 

In summary, this paper demonstrates 
that a reliable methodology can be used to 
measure Young's modulus and step height. 
Furthermore, a micro and nano technology 
(MNT) 5-in-l standard reference material 
(SRM) can be used by industry to compare 
their in-house measurements using this 
methodology with NIST measurements 
thereby validating their use of the 
documentary standards. 
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List of Terms 

cantilever 

fixed-fixed beam 



a test structure that consists of a 
freestanding beam that is fixed at 
one end ' 

a test structure that consists of a 
freestanding beam that is fixed at 
both ends ' 



in-plane length (or deflection) measurement 

the experimental determination of 
the straight-line distance between 
two transitional edges in a MEMS 
device ' 

interferometer a non-contact optical instrument 

used to obtain topographical 3-D 
data sets ' 



Reprinted, with permission, from ASTM E 2444 Terminology Relating to Measurements Taken on Thin, Reflecting Films, copyright ASTM 
International, 100 Barr Harbor Drive, West Conshohocken, PA 19428. 
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residual strain 



residual stress 



step height 



(residual) strain 
gradient 



(residual) stress 
gradient 



test structure 



thickness 

vibrometer 
Young's modulus 



in a MEMS process, the amount of 
deformation (or displacement) per 
unit length constrained within the 
structural layer of interest after 
fabrication yet before the con- 
straint of the sacrificial layer (or 
substrate) is removed (in whole or 
in part) ' 

the remaining forces per unit area 
within the structural layer of 
interest after the original cause(s) 
during fabrication have been 
removed yet before the constraint 
of the sacrificial layer (or sub- 
strate) is removed (in whole or in 
part) 2 

the distance in the z-direction that 
an initial, flat, processed surface 
(or platform) is to a final, flat, 
processed surface (or platform) 2 
a through-thickness variation (of 
the residual strain) in the structur- 
al layer of interest before it is 
released ' 

a through-thickness variation (of 
the residual stress) in the structur- 
al layer of interest before it is 
released 2 

a component (such as, a fixed- 
fixed beam or cantilever) that is 
used to extract information (such 
as, the residual strain or the strain 
gradient of a layer) about a fabri- 
cation process ' 

the height in the z-direction of one 
or more designated thin-film 
layers 2 

an instrument for non-contact 
measurements of surface motion 2 
a parameter indicative of material 
stiffness that is equal to the stress 
divided by the strain when the 
material is loaded in uniaxial 
tension, assuming the strain is 
small enough such that it does not 
irreversibly deform the material 2 



Reprinted, with permission, from SEMI MS2, MS3, and MS4 
copyright Semiconductor Equipment and Materials International, 
Inc. (SEMI)© 2010, 3081 Zanker Road, San Jose, CA 95134, 
www.semi.org. 



List of Symbols 

For Young's modulus measurements and calculations: 

ft = viscosity of the ambient surround- 

ing the cantilever. [SEMI MS4] 

p = density of the thin film layer. [SEMI MS4] 

E = calculated Young's modulus value of the 

thin film layer. [SEMI MS4] 

^damped = calculated Young's modulus value ob- 

tained from the average resonance fre- 
quency of a fixed-fixed beam assuming 
clamped-clamped boundary conditions. 
[SEMI MS4] 

E- mit = initial estimate for the Young's modulus 

value of the thin film layer. [SEMI MS4] 

^simple = calculated Young's modulus value ob- 

tained from the average resonance fre- 
quency of a fixed-fixed beam assuming 
simply supported boundary conditions 
at both supports. [SEMI MS4] 

f cm = average undamped resonance frequen- 

cy of the cantilever. [SEMI MS4] 

f m = average resonance frequency of the 

fixed-fixed beam. [SEMI MS4] 

L cm = suspended cantilever length. [SEMI MS4] 

Zfl, = suspended fixed-fixed beam length. 

[SEMI MS4] 

Q = oscillatory quality factor of the 

cantilever. [SEMI MS4] 

t = thickness of the thin film layer. [SEMI 

MS4] 

W c<m = suspended cantilever width. [SEMI MS4] 

For combined standard uncertainty calculations of 
Young's modulus measurements: 

O^ = one sigma uncertainty of the value of fi. 

[SEMI MS4] 
G p = one sigma uncertainty of the value of p. 

[SEMI MS4] 
<7 freq = one sigma uncertainty of the value of 

f cm . [SEMIMS4] 
a L = one sigma uncertainty of the value of 

Z can . [SEMIMS4] 
°fhick = one sigma uncertainty of the value of t. 

[SEMI MS4] 
a w = one sigma uncertainty of the value of 

W cm . [SEMI MS4] 
E misyi = maximum Young's modulus value as 

determined in an uncertainty calcula- 
tion. [SEMI MS4] 
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^damp 



' livq 



E miR = minimum Young's modulus value as de- 

termined in an uncertainty calculation. 
[SEMI MS4] 

f msol = frequency resolution for the given set of 

measurement conditions. [SEMI MS4] 

p m = estimated percent difference between 

the damped and undamped resonance 
frequency of the cantilever. [SEMI MS4] 

u p = component in the combined standard 

uncertainty calculation for Young's 
modulus that is due to the uncertainty 
of p. [SEMIMS4] 

u c = combined standard uncertainty value 

(that is, the estimated standard deviation 
of the result). [SEMIMS4] 

u cE = the combined standard uncertainty for a 

Young's modulus measurement. 
= component in the combined standard 
uncertainty calculation for Young's 
modulus that is due to damping. [SEMI 
MS4] 
= the standard uncertainty for a Young's 
modulus measurement as obtained from 
a fixed-fixed beam. 
= component in the combined standard 
uncertainty calculation for Young's 
modulus that is due to the measurement 
uncertainty of f can . [SEMI MS4] 

Ur = component in the combined standard 

uncertainty calculation for Young's 
modulus that is due to/ resol . [SEMI MS4] 

u L = component in the combined standard 

uncertainty calculation for Young's 
modulus that is due to the measurement 
uncertainty of Z can . [SEMIMS4] 

"thick = component in the combined standard 

uncertainty calculation for Young's 
modulus that is due to the measurement 
uncertainty of t. [SEMI MS4] 

For calibration of step height measurements: 

a cert = the one sigma uncertainty of the cali- 

brated physical step height used for 
calibration. 

cert = the certified value of the double-sided 

physical step height used for calibration. 
[SEMI MS2] 

z drift = the calibrated positive difference 

between the average of the six calibra- 
tion measurements taken before the data 
session (at the same location on the 



-repeal 



physical step height used for calibra- 
tion) and the average of the six calibra- 
tion measurements taken after the data 
session (at the same location on the 
physical step height used for calibra- 
tion). [SEMI MS2] 
= over the instrument's total scan range, 
the maximum percent deviation from 
linearity, as quoted by the instrument 
manufacturer (typically less than 3 %). 
[SEMI MS2] 
= the maximum of two calibrated values; 
one of which is the positive calibrated 
difference between the minimum and 
maximum values of the six calibration 
measurements taken before the data 
session (at the same location on the 
physical step height used for calibra- 
tion) and the other is the positive 
calibrated difference between the 
minimum and maximum values of the 
six calibration measurements taken after 
the data session (at the same location on 
the physical step height used for calibra- 
tion). [SEMI MS2] 

z = calibrated average of the twelve calibra- 

tion measurements taken before and 
after the data session at the same loca- 
tion on the physical step height used for 
calibration. [SEMI MS2] 

z 6 = the calibrated average of the six calibra- 

tion measurements from which z repeat is 
found. [SEMI MS2] 

For step height measurements and calculations: 



J platNXt 



-'roughAiY 



platNrD 



the calibrated standard deviation of the 
data from trace "f on platNX. [SEMI 
MS2] 

the calibrated surface roughness of 
platNX measured as the smallest of all 
the values obtained for G platNXt ; however, 
if the surfaces of the platforms (includ- 
ing the reference platform) all have 
identical compositions, then it is meas- 
ured as the smallest of all the standard 
deviation values obtained from data 
traces "a," "b," and "c" along these 
platforms. [SEMI MS2] 
the average of the calibrated reference 
platform height measurements taken 
from multiple data traces on one step 
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height test structure, where N is the test 
structure number ("1," "2," "3," etc.), 
r indicates it is from a reference plat- 
form, and D directionally indicates 
which reference platform (using the 
compass indicators "N," "S," "E," or 
"W" where "N" refers to the reference 
platform designed closest to the top of 
the chip). [SEMI MS2] 

platNrDt = a calibrated reference platform height 
measurement from one data trace, where 
TV is the test structure number ("1," "2," 
"3," etc.), r indicates it is from a refer- 
ence platform, D directionally indicates 
which reference platform (using the 
compass indicators "N," "S," "E," or 
"W" where "N" refers to the reference 
platform designed closest to the top of 
the chip), and t is the data trace ("a," 
"b," "c," etc.) being examined. [SEMI 
MS2] 

platNX = the calibrated platform height measure- 
ment, where N is the test structure 
number ("1," "2," "3," etc.) and X is the 
capital letter (or "r" is used if it is the 
reference platform) associated with the 
platform ("A," "B," "C," etc.) as 
lettered starting with "A" for the 
platform closest to platNrW or platNrS. 
[SEMI MS2] 

platNXt = a calibrated platform height measure- 
ment from one data trace, where N is the 
test structure number ("1," "2," "3," 
etc.), X is the capital letter associated 
with the platform ("A," "B," "C," etc.) 
as lettered starting with 'A" for the 
platform closest to platNrW or platNrS, 
and t is the data trace ("a," "b," "c," etc.) 
being examined. [SEMI MS2] 

stepN x M Y = the calibrated step height measurement 
taken from two different step height test 
structures (N and M) on the same test 
chip, which is equal to the final platform 
height minus the initial platform height, 
where the step is from the initial plat- 
form to the final platform and where X is 
the capital letter associated with the 
initial platform from test structure 
number TV, and Y is the capital letter 
associated with the final platform from 
test structure number M. [SEMI MS2] 

stepNxY = the average of the calibrated step height 



measurements taken from multiple 
data traces on one step height test struc- 
ture, where N is the number associated 
with the test structure, X is the capital 
letter associated with the initial plat- 
form (or "r" is used if it is the reference 
platform), Y is the capital letter associat- 
ed with the final platform (or "r" is used 
if it is the reference platform), and the 
step is from the initial platform to the 
final platform. [SEMI MS2] 
stepNxY, = a calibrated step height measurement 
from one data trace on one step height 
test structure, where N is the number 
associated with the test structure, X is 
the capital letter associated with the 
initial platform (or "r" is used if it is the 
reference platform), Y is the capital 
letter associated with the final platform 
(or "r" is used if it is the reference plat- 
form), t is the data trace ("a," "b," "c," 
etc.) being examined, and the step is 
from the initial platform to the final 
platform. [SEMI MS2] 

For combined standard uncertainty calculations of step 
height measurements: 

u c = the combined standard uncertainty value 

(i.e., the estimated standard deviation of 
the result). [SEMI MS2] 

u cert = the component in the combined standard 

uncertainty calculation that is due to the 
uncertainty of the value of the physical 
step height used for calibration. [SEMI 
MS2] 

w cSH = the combined standard uncertainty or a 

step height measurement. 

u Mft = the uncertainty of a measurement due to 

the amount of drift during the data 
session. 

"linear = the uncertainty of a measurement due to 

the deviation from linearity of the data 
scan. [SEMI MS2] 

u l p mnx = the component in the combined standard 

uncertainty calculation for platform 
height measurements that is due to the 
measurement uncertainty across the 
length of platNX, where the length is 
measured perpendicular to the edge of 
the step. [SEMI MS2] 

u Lstep = the component in the combined standard 
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' repeal 



l WplatNX 



uncertainty calculation for step height 
measurements that is due to the 
measurement uncertainty of the step 
height across the length of the step, 
where the length is measured perpendi- 
cular to the edge of the step. [SEMI 
MS2] 
'-tphtNx = the component in the combined standard 

uncertainty calculation or step height 
measurements obtained from two step 
height test structures that is due to the 
uncertainty of the platform height 
measurement for platNX. [SEMI MS2] 

= the uncertainty of a measurement due to 
the repeatability of a measurement. 
[SEMI MS2] 

= the component in the combined standard 
uncertainty calculation for platform 
height measurements that is due to the 
measurement uncertainty across the 
width of platNX, where the width is 
measured parallel to the edge of the 
step. [SEMI MS2] 
u wste P = the component in the combined standard 

uncertainty calculation for step height 
measurements that is due to the meas- 
urement uncertainty of the step height 
across the width of the step, where the 
width is measured parallel to the edge of 
the step. [SEMI MS2] 

For residual stress and stress gradient calculations: 

e r = residual strain of the thin film layer. 

[SEMI MS4] 

<7 g = stress gradient of the thin film layer. 

[SEMI MS4] 

°gmax = maximum stress gradient value as deter- 

mined in an uncertainty calculation. 
[SEMI MS4] 

(Tgnun = minimum stress gradient value as deter- 

mined in an uncertainty calculation. 
[SEMI MS4] 



a r 

ffrm 



1 t:r( Gr) 



*E(<jg) 



l aar) 



II 



..KKTL') 



residual stress of the thin film layer. 
[SEMI MS4] 

maximum residual stress value as deter- 
mined in an uncertainty calculation. 
[SEMI MS4] 

minimum residual stress value as deter- 
mined in an uncertainty calculation. 
[SEMI MS4] 

strain gradient of the thin film layer. 
[SEMI MS4] 

component in the combined standard 
uncertainty calculation for residual 
stress that is due to the measurement 
uncertainty of £ t . [SEMI MS4] 
combined standard uncertainty value for 
residual strain. 

combined standard uncertainty value for 
stress gradient. 

combined standard uncertainty value for 
residual stress. 

combined standard uncertainty value for 
strain gradient. 

component in the combined standard 
uncertainty calculation for stress gradi- 
ent that is due to the measurement 
uncertainty of E. [SEMI MS4] 
component in the combined standard 
uncertainty calculation for residual 
stress that is due to the measurement 
uncertainty of E. [SEMI MS4] 
component in the combined standard 
uncertainty calculation for stress gradi- 
ent that is due to the measurement 
uncertainty of s g . [SEMI MS4] 



For round robin measurements: 



the number of reproducibility or repea- 
tability measurements, 
the average combined standard uncer- 
tainty value for the reproducibility or 
repeatability measurements. It is equal 
to the sum of the u c values divided by n. 
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1. Introduction 

Microelectromechanical systems (MEMS) foundaries 
are emerging from the economic downturn in a strong 
position [1]. Applications for MEMS demand robust 
reliability [2], which necessitates the consideration of 
reliability at the earliest stages of product development 
[3]. The U.S. Measurement System (USMS) stated that 
there are insufficient measurement data and measure- 
ment methods to generate that data to adequately 
characterize a MEMS device and its dependence on 
fabrication processes [4]. In response to this, SEMI has 
published two MEMS standard test methods; one for 
measuring Young's modulus [5] and one for measuring 
step heights [6]. This paper presents the technical basis 
for the standards and the data from a MEMS Young's 
modulus and step height round robin experiment, 
which was used to validate these standards. These 
standards are expected to facilitate commerce in 
MEMS technologies and improve manufacturing yields 
by decreasing interlaboratory differences in these 
measurements. 

The first MEMS standards were approved by ASTM 
in 2002 then validated in 2005 [7-9] with round robin 
precision and bias data. These ASTM test methods are 
used for measuring in-plane lengths or deflections, 
residual strain, and strain gradient. They employ a non- 
contact measurement approach using an optical inter- 
ferometer. 3 

The SEMI test method for measuring Young's 
modulus that is reported in this paper uses an optical 
vibrometer, stroboscopic interferometer, or comparable 
instrument, and the SEMI test method for step height 
measurements also reported in this paper uses an 
optical interferometer or comparable instrument. 
Therefore, a stroboscopic interferometer can be used 
for all five standard test methods for MEMS (the three 
ASTM test methods and the two SEMI test methods). 

Test structures for use with all five standard test 
methods are incorporated on the MEMS Young's 
Modulus and Step Height Round Robin Test Chip, 
which is shown in Fig. 1 . This chip is also the prototype 
for a micro and nano technology (MNT) 5-in-l 
standard reference material (SRM) [10]. Using the 
MNT 5-in-l SRM, companies will be able to compare 
their in-house measurements taken on the SRM with 
NIST measurements thereby validating their use of the 
documentary standards. 



In this paper, commercial equipment, instruments, or processes 
may be identified. This does not imply recommendation or endorse- 
ment by NIST, nor does it imply that the equipment, instruments, or 
processes are the best available for the purpose. 



The first test method (SEMI MS4) [5] to be present- 
ed in this paper is on Young's modulus measurements 
of thin, reflecting films. As presented in [11], there are 
many Young's modulus methods with values that have 
been reported in the literature. The method reported 
here is based on the textbook treatment [12] of the aver- 
age resonance frequency of a single-layered cantilever 
oscillating out-of-plane (recommended and shown in 
Fig. 2) or based on the average resonance frequency of 
a single layered fixed-fixed beam oscillating out-of- 
plane. The resonance frequency method was chosen 
due to a) its wide acceptance in the community, b) the 
involvement of simple, non-contact, non-destructive 
measurements, and c) its ability to obtain measure- 
ments from test structures on a single test chip along- 
side other test structures. A piezoelectric transducer 
(PZT) is used to create the out-of-plane excitations; 
however, a PZT is not the only means available for 
excitation (for example, thermal excitation [13-15] is 
possible). Given the value for Young's modulus, 
residual stress and stress gradient calculations are 
possible. 

Young's modulus measurements are an aid in the 
design and fabrication of MEMS devices and inte- 
grated circuits (ICs). For example, high values of 
residual stress can lead to failure mechanisms in 
ICs such as electromigration, stress migration, and 
delamination. So, methods for its characterization 
are of interest for IC process development and monitor- 
ing in order to improve the yield in complementary 
metal oxide semiconductor (CMOS) fabrication 
processes [16]. 

The second test method (SEMI MS2) [6] to be pre- 
sented in this paper is for step height measurements of 
thin films. A step height test structure, such as shown in 
Fig. 3(a) with its cross section given in Fig. 3(b), is 
used for these measurements. Multiple 2-D data traces 
along the top of this test structure, such as traces "a," 
"b," and "c" in Fig. 3(a), are obtained. A sample 2-D 
data trace is given in Fig. 3(c). Data averages and 
standard deviations from the 2-D data traces along the 
pertinent platforms are recorded. 

For data obtained from one step height test structure, 
the difference in the platform heights involved in the step 
for each 2-D data trace is calculated. The step height, 
stepN^, [or stepl AB as shown in Fig. 3(b)] is the average 
of these values from the different 2-D data traces. 

Step height measurements can also be obtained from 
two step height test structures. The calculations are 
slightly different, and this approach is not recommend- 
ed due to higher resulting combined standard 
uncertainty values [17]. 
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Fig. 1. The MEMS Young's Modulus and Step Height Round Robin Test Chip. 



dimensional 
marker 




' etch stop 



H 



exposed silicon 
for etch 
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dimensional 
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etch stop 



Fig. 2. For a cantilever test structure on the test chip shown in Fig. 1 (a) a design rendition and (b) a cross section. In Fig. 2(a), the central 
structure (silicon dioxide and silicon nitride atop the silicon substrate) will become an all oxide cantilever after the silicon nitride is removed and 
the exposed raw silicon on three of the structures sides and beneath it is etched, as shown in Fig. 2(b). The dimensional markers are typically made 
of polysilicon or metal encapsulated in oxide and can be used to obtain a more accurate measurement of the cantilever length after the 
post-processing etch. 
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Fig. 3. For the step height test structure, on the test chip shown in Fig. 1, with a 0° 
orientation (a) a design rendition showing 2-D data traces "a," "b," and "c," (b) its cross 
section, (where aa indicates active area, fox indicates field oxide, pi indicates polyl, m2 
indicates metal2, and gl indicates glass), and (c) a sample 2-D data trace. 



Step height measurements can be used to determine 
thin film thickness values [18,19]. Thickness measure- 
ments are an aid in the design and fabrication of MEMS 
devices and can be used to obtain thin film materials 
parameters, such as Young's modulus [5,16]. 

The SEMI test methods for Young's modulus and step 
height measurements were used in the MEMS Young's 
modulus and step height round robin experiment to 
determine the repeatability of a measurement as well as 
to see if independent laboratories could reproduce these 
measurements without introducing a bias. 



For the round robin, test chips were fabricated on the 
same processing run. One test chip (the design of which 
is shown in Fig. 1) was delivered to each laboratory, 
and measurements were taken using the procedures in 
the SEMI test methods. The data from these measure- 
ments were analyzed using a NIST web-based program 
[10] that also verifies the data. Five laboratories 
(including NIST) participated in the round robin, 
resulting in eight independent data sets (due to different 
equipment setups) for Young's modulus and seven 
independent data sets (due to different equipment 
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setups) for step heights. The results from this round 
robin are reported in this paper. 

Section 2 presents the packaged Round Robin 
Test Chip. Section 3 presents the technical basis for 
the Young's modulus and step height standard test 
methods. Section 4 presents the uncertainty equations. 
And, Sec. 5 presents the results of the MEMS Young's 
modulus and step height round robin experiment 
followed by the conclusions in Sec. 6. Instrument 
specifications are given in Appendix A (Sec. 7). 
Supplemental material for Young's modulus measure- 
ments can be found in Appendix B (Sec. 8), and supple- 
mental material for step height measurements can be 
found in Appendix C (Sec. 9). 

2. Packaged Round Robin Test Chip 

The MEMS Round Robin Test Chip was fabricated 
on a 1.5 urn CMOS process available through MOSIS 
[20]. The design for this test chip is depicted in Fig. 1. 
The design file (in GDS-II format) for this chip can be 
downloaded from the NIST Semiconductor Electronics 
Division (SED) MEMS Calculator website [10]. 

For the round robin chip design shown in Fig. 1 , in a 
number of places one mechanical layer is fabricated 
into suspended structures such as cantilevers and 
fixed-fixed beams. This layer consists of all oxide; 
namely, the field oxide, the deposited oxide before and 
after the first metal (ml) deposition, and the glass (gl) 
layer. [The nitride cap (present atop the glass layer 
when the chips are received from MOSIS) was 
removed after fabrication using a CF 4 + 2 etch before 
a post-processing XeF 2 etch that released the beams.] 

As seen in Fig. 1, the test chip contains six groupings 
of test structures with the following labels: 

1 . Young's Modulus, 

2. Residual Strain, 

3 . Strain Gradient, 

4. Step Height, 

5. In-Plane Length, and 

6. Certification Plus. 

However, for the MEMS Young's modulus and 
step height round robin experiment, we were only con- 
cerned with the first and fourth groupings of test struc- 
tures, for Young's modulus and step height measure- 
ments, respectively. 

In the first grouping of test structures on the Round 
Robin Test Chip shown in Fig. 1, Young's modulus 
measurements were made. Cantilever and fixed-fixed 



beam test structures were provided for this purpose with 
25 cantilevers grouped above 25 fixed-fixed beams on 
this chip; however, for the round robin we were only 
concerned with the cantilevers, such as that shown in 
Fig. 2. Configurations for the cantilevers on the chip 
shown in Fig. 1 are given in Table 1. (See Sec. 3.1 for the 
rationale for the chosen cantilever dimensions.) 

Table 1. Cantilever configurations for Young's modulus measure- 
ments 



Test Structure 
Width (in um) 
Length (in um) 
Orientation 
Mechanical Layer 
Quantity of Beams 



Cantilevers 

28 

200, 248, 300, 348, and 400 

0° 

Oxide 

5 of each length (or 25 beams) 



As shown in Table 1, the cantilever design lengths 
are 200 um, 248 um, 300 um, 348 um, and 400 um. 
The length of the cantilever (in micrometers) is given at 
the top of each column of cantilevers in Fig. 1 follow- 
ing the column number (i.e., 1 to 5). The beams are 
designed at only a 0° orientation. 4 There are five 
cantilevers designed at each length. Therefore, there are 
25 oxide cantilevers with a 0° orientation. 

In the fourth grouping of test structures on the Round 
Robin Test Chip shown in Fig. 1 , step height measure- 
ments are made. Step height test structures, such as 
shown in Fig. 3, are provided for this purpose for a 
metal2-over-polyl (m2-over-pl) step from active area 
(aa) to field oxide (fox). The surrounding reference 
platform consists of the deposited oxides and metal2 
over active area. The metal2 thickness is approximately 
1.0 um. 

There are four orientations (0°, 90°, 180°, and 270°) 
of the step height test structure shown in Fig. 3(a), and 
these orientations are grouped in quads. For the Round 
Robin Test Chip in Fig. 1 there are three quads, one of 
which is shown in Fig. 4. The quad number is given in 
the center of the quad. Therefore, quad "2" (out of 3) is 
shown in Fig. 4. 

After the chips were received from MOSIS, they 
were post-fabricated in a class 100 clean room (to 
remove the nitride cap with a CF 4 + 2 etch and to 



A 0° orientation implies that the length of the beam is parallel to 
the x-axis of the test chip, the axes of which are shown in Fig. 1 and 
again in Fig. 2, with the connection point of the cantilever having 
a smaller x-value than the x-values associated with the suspended 
portion of the cantilever. 
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AA-FO 




V 



°-*.V 



FO-AA 



Fig. 4. A design rendition of quad "2" from the step height grouping 
of test structures in Fig. 1. 



release the beams with a XeF 2 etch), packaged in a 
laboratory environment, then stored in a dust-free N 2 
atmosphere. 

Each round robin participant received a packaged 
Round Robin Test Chip, as shown in Fig. 5. The pack- 
aged part was put together in the following way: 

1. Starting with a hybrid package with a pin 
arrangement similar to that shown in Fig. 5, the 
PZT was secured to the top of the chip cavity 
using two thin layers of low stress, non-conduct- 



ing epoxy. (The first layer of epoxy ensures that 
there will not be a conducting path between the 
package and the PZT. ) 

2. The PZT has the following properties: 

a. The dimensions of the PZT are approxi- 
mately 5 mm by 5 mm and 2 mm in height, 

b. It is provided with a red and a black wire, 

c. It can achieve a 2.2 um (± 20 %) displace- 
ment at 100 V, 

d. It has an electrical capacitance of 250 nF 
( ± 20 %), and 

e. It has a resonance frequency greater than 
300 kHz, at which or above which it shall 
not be operated because it could damage the 
PZT. 

3. The two PZT wires were secured to their respec- 
tive package connections. 

4. The Round Robin Test Chip was secured to the 
top of the PZT using two thin layers of a low 
stress non-conducting epoxy. (The first layer of 
epoxy ensures that there will not be a conducting 
path between the PZT and the test chip. ) 

5. The lid (or can) was placed on top of the pack- 
age to protect the chip, and the can was secured 
to the package with tape before shipment. 



to vibro meter or 

comparable 
instrument 




non-conductive 
epoxy 



Fig. 5. The packaged Round Robin Test Chip. 
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To take measurements on the Round Robin Test 
Chip, the can is carefully removed. Young's modulus 
and step height measurements can now be taken. For 
Young's modulus measurements, to operate the PZT, 
the red wire is driven with a voltage that is positive rel- 
ative to the black wire. To ensure that you have suc- 
cessfully connected to the PZT, when activated at 10 V 
and 7000 Hz, the resulting PZT vibration is barely 
audible. 



If these are damped measurements they are called 
/dampedi, /damped and /ctapeda, respectively. For each 
damped frequency (4^ „ / damped2 , andf &amptdi ), a corre- 



sponding undamped frequency (f u 



b inn 



undamped!' J Lmdamped2' 



and 

Andam P ed3. respectively) is calculated using the equation 
below: 



/;, 



/d 



dumped/, 



undamped n 



1 
4cf 



(2) 



3. Young's Modulus and Step Height 
Measurements 

The technical basis for the standard test method on 
Young's modulus measurements is given in Sec. 3.1, 
and the technical basis for the standard test method on 
step height measurements is given in Sec. 3.2. 



where the n in the subscript of^™^ and/ u 



/ undamped?; 



is 1, 



2, or 3 and where Q is the oscillatory quality factor of 
the cantilever as given by the following equation [21]: 



Q- 



24\i 



( t 


\ 2 


L 



(3) 



3. 1. Young's Modulus Measurements 

The Young's modulus of a single layer is obtained 
from resonance frequency measurements of a can- 
tilever comprised of that layer (such as shown in Fig. 2) 
or from resonance frequency measurements of a fixed- 
fixed beam comprised of that layer. To determine an 
estimate for the fundamental resonance frequency of a 
cantilever, / caninit , the following equation (as derived in 
Sec. 8.1) is used: 5 



J ca 



E..J 2 



38.330p4, 



(1) 



where W ca „ is the suspended cantilever width and jx is 
the viscosity (in air, yt = 1.84 X 10 s Ns/m 2 at 20 °C). 

The cantilever dimensions for the Round Robin Test 
Chip in Fig. 1, as specified in Table 1, were chosen 
such that 5 urn < W can < 40 urn, W cm > t, and L can » t 
where t = 2.743 um, as determined by the electro- 
physical technique [18]. Data Sheet Tl [10] can be 
used to calculate t. In addition, the cantilever dimen- 
sions were chosen to achieve a) an estimated reso- 
nance frequency between 10 kHz and 75 kHz using 
Eq. (1) with the assumptions that is,„„=70GPa and 
p = 2.2 g/cm 3 , b) a Q value above 30 using Eq. (3), and 
c) a value less than 2 % for/> rf# as given by the follow- 
ing equation: 



where E Mt is the initial estimate for the Young's modu- 
lus value of the thin film layer, t is the thickness, p is 
the density, and L cm is the suspended length of the 
cantilever. 

Measurements are taken at frequencies which 
encompass f cmMt , and an excitation-magnitude versus 
frequency plot is obtained from which the resonance 
frequency is found. For a given cantilever, three meas- 
urements of resonance frequency are obtained (namely, 
/measi./me^, and/ meas3) . If these are undamped measure- 
ments (e.g., if they are performed in a vacuum) they 
are called Andampedb /undamped2> an d /undamped3> respectively. 



By inserting the inputs into the correct locations on the appropri- 
ate NIST Web page (http://www.eeel.nist.gov/812/test-structures/ 
MEMSCalculator.htm), the calculations can be performed on-line in 
a matter of seconds. 



.1' 

J u: 



damped n 



ndamped n 



100° 



1-1 



1 



\Q l 



100%. (4) 



See Table 2 for the calculations of f ca „ Mt , Q, and^ d ,yfor 
the chosen dimensions. 

Table 2. Calculations of / cammt , Q, and p di tr for the designed 
cantilever lengths ( W cm = 28 um) 



^can 

(um) 


J caninit 
(kHz) 


Q 


Pdiff 

(%) 


200 


62.5 


148.0 


0.0006 


248 a 


40.6 


96.3 


0.0013 


300 


27.8 


65.8 


0.0029 


348" 


20.6 


48.9 


0.0052 


400 


15.6 


37.0 


0.0091 



These values were chosen in order to design on a 0.8 um grid to 
simplify the interface with MOSIS and the fabrication facility. 
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Also, to ensure that the resonance frequency of the 
cantilever is not altered by squeeze film or other damp- 
ing phenomena, 6 the cantilever should be suspended 
high enough above the underlying layer such that its 
motion is not altered by the underlying layer. 7 In other 
words, the gap, d, between the suspended cantilever 
and the underlying layer should adhere to the following 
equation [22]: 



d> 



W„ 



(5) 



Therefore, if W can = 28 urn, the gap between the 
suspended cantilever and the underlying layer should 
be at least 9.3 um. 

The average undamped resonance frequency, f can , is 
calculated from the three undamped resonance frequen- 
cies using the following equation: 



.L 



J undamped 1 ./undamped2 Aindamped3 



(6) 



Given the measured value fovf can , the Young's modulus 
value, E, is calculated as follows: 



A%n 2 pflLl 



38.330p/" 2 £ 

r Jean cat 



1.875 4 



(7) 



which assumes clamped-free boundary conditions and 
no undercutting of the beam. The derivation of this 
equation is presented in Sec. 8.1. The combined 
standard uncertainty for E, or u cE , is given in Sec. 4.2. 
Residual stress and stress gradient equations for this 
thin film layer can be found in Sec. 8.2. Consult 
Sec. 8.3 for Young's modulus measurements obtained 
from fixed-fixed beams. 



3.2 Step Height Measurements 

Step height measurements can be taken from either 
one or two step height test structures; however, 
measurements from one step height test structure are 



recommended due to lower resulting values for the 
combined standard uncertainty, w cSH [17]. 

If one step height test structure is used to obtain a 
step height measurement, three 2-D data traces ["a," 
"b," and "c," as shown in Fig. 3(a)] are taken along the 
top of the test structure, a cross section of which is 
given in Fig. 3(b). A sample 2-D data trace is shown in 
Fig. 3(c). All height measurements are with respect to 
the height of the surrounding reference platform that is 
used to level and zero the data. For generic test struc- 
ture "N" with platforms labelled "A"' and "Y," the indi- 
vidual platform height measurements (namely, 
platNXa, platNXb, platNXc, platNYa, platNYb, and 
platNYc) and the standard deviations {o plamx „ <J platNXb , 



GplatNXo GplatNYzi GplatNYb, atl d ^ptaATc) from the tWO plat- 

forms involved in the step in data traces "a," "b," and 
"c" are recorded. 8 If the test structure in Fig. 3(a) is 
called test structure "1," then for the step in test struc- 
ture "1" from platform "A" to platform "B," the plat- 
form height measurements would be platlAa, platlAb, 
plat\ Ac, platlBa, platlBb, and platlBc, and the 
standard deviations would be a planM , a plallAb , a platlAc , 
<Vtma, entomb, and <7 ptolBc . Therefore, 12 measurements 
are obtained (6 from platform "A" and 6 from platform 
"B"). 

For each 2-D data trace, the difference in the plat- 
form heights involved in a step (in general, stepN^ 9 
is calculated using the following equation: 



stepN^ = platNYt - platNXt 



(8) 



where t is the data trace ("a," "b," "c," etc.) being 
examined. For the step shown in Fig. 3(b) from 
platform 'A" to platform "B," the equations are: 



and 



stepl^^ = platlBa — platlAa , 
stepl^ = platlBb — platl Ab , 

step\ ABc = platlBc - platl Ac . 



(9) 
(10) 

(11) 



Squeeze film and other damping phenomena routinely lead to amplitude dependent resonance frequencies and shifts in the natural frequency 
of the system, which may limit the accuracy of the technique. 

The damping may not be present in bulk-micromachining processes because it is dependent upon the depth of the cavity and the vicinity of the 
sides of the cavity to the beam. (The damping phenomena are expected to be present in surface micromachining processes without the use of a 
backside etch, unless the measurement is performed in a vacuum.) 

Consult the List of Symbols for the nomenclature used for platNXt and O p i atNX t ■ 

Consult the List of Symbols. 
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The step height, stepNxy, is the average of the 
values from the different 2-D data traces as given 
below: 



stepN x 



stepN^ + stepN^ + stepN ^ ^ 



or for the step shown in Fig. 3(b), the step height, 
stepl m , is: 



stepl A 



^P 1 ^ + Ste P 1 ABb +Ste P 1 ABc / 13) 



The combined standard uncertainty, m cSH , for stepNxy 
is given in Sec. 4.3. Consult Sec. 9. 1 for a step height 
measurement obtained from two step height test 
structures. 



4. Uncertainty Equations 

In this section, the equations used to determine the 
values of the combined standard uncertainty [17], u c , 
are presented. Sec. 4.1 presents the basic combined 
standard uncertainty equation, and Sec. 4.2 and Sec 4.3 
present the more specific uncertainty equations for 
Young's modulus and step height measurements, 
respectively. 

4.1 Combined Standard Uncertainty Equation 

The combined standard uncertainty, u c , [17] is calcu- 
lated as the estimated standard deviation of the result. It 
is equal to the square root of the sum of the squares of 
the uncertainty components where each component 
must have the same units as u c (e.g., GPa). 

For the case of three sources of uncertainty, the 
uncertainty equation would be as follows: 



4.2 Young's Modulus Uncertainty Equations 

In this section, a combined standard uncertainty 
equation is presented for use with Young's modulus 
measurements (u cE ) as obtained from resonance fre- 
quency measurements from a cantilever. For these 
measurements, six sources of uncertainty are identified 
with all other sources considered negligible. The six 
sources of uncertainty are the uncertainty of the thick- 
ness (u tUck ), the uncertainty of the density (u p ), the 
uncertainty of the cantilever length (u L ), the uncertain- 
ty of the average resonance frequency (u^), the uncer- 
tainty due to the frequency resolution (w^ o/ ), and the 
uncertainty due to damping (u damp ). As such, the com- 
bined standard uncertainty equation for u cE (as calculat- 
ed in SEMI Test Method MS4 [5] and in Data Analysis 
Sheet YM. 1 [10]) with six sources of uncertainty is as 
follows: 



: V^ 



2 
thick 



+ u p + u L + u) nq + u fresol + ul mp 



(15) 



In determining the combined standard uncertainty, a 
Type B evaluation [17] (i.e., one that uses means other 
than the statistical Type A analysis) is used for each 
source of error. Table 3 gives sample values for each of 
these uncertainty components. 

The uncertainty for u lhjck is determined from calculat- 
ed minimum and maximum Young's modulus values 
(namely, E min and E miD . , respectively). It is derived using 
the extremes of values expected for the cantilever 
thickness as given below: 



38.330p/j„r 



2 T 4 

'am 



and 



+ 3<w) 



38.330p/lLl 



(16) 



(17) 



yjU, 



(14) 



where u x is the uncertainty component due to the first 
source of uncertainty, u 2 is the uncertainty component 
due to the second source of uncertainty, and u 3 is due to 
the third source of uncertainty. Additional terms may be 
added under the square root sign for additional sources 
of uncertainty. 



where O thick is the one sigma uncertainty of the value of 
t, which is found using the electro-physical technique 
[18]. Data Sheet Tl [10] can be used to calculate (J thick . 
With 99.7 % confidence, assuming a Gaussian distribu- 



tion (and assuming u 



p, U L , Ufreq, 



Ufr eso ,, and u damp equal 



zero), the value for E lies between E n 
Therefore, u Mclc is calculated as follows: 



and E„ 



E - E 

max min 



(18) 
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Table 3. Sample Young's modulus uncertainty values (assuming E jnit = 70 GPa) 

uncertainty 

"thick 

(using t = 2.1 A3 um and CT^^ = 0.058 urn) 

Up 

(using p = 2.2 g/cm 3 and G„ = 0.05 g/cm 3 ) 
(using Z cim = 300 um and a L = 0.2 um) 



"L 



"freq 



"diimp 



"cE 



3» C £ 



(using/ measl = 26.82625 kHz, 
/ meas2 = 26.8351 kHz, 
and / meas3 = 26.8251 kHz) 



(using / resol = 1.25 Hz) 



Type A or Type B 


value (in GPa) 


Type B 


2.8 


Type B 


1.5 


Type B 


0.17 


Type B 


0.027 



(using W cm = 28 um and O w = 0.1 um 

and using n = 1.84 X 10~ 5 Ns/m 2 and a fl = 0.01 X 10~ 5 Ns/m 2 ) 



V"ick + "p + u l + u L q + «Li + » 



damp 



Type B 
Type B 



0.0018 



0.0004 



3.2 



9.5 



This 3w c £ uncertainty is plotted in Fig. 6 with the repeatability data point corresponding to the first cantilever with length of 300 |j,m. 



In the same way, to determine the uncertainty com- 
ponent, u p , the calculated minimum and maximum 
Young's modulus values (namely, is min and E mia , respec- 
tively) are given below: 

38.330 (p-3a p )/ c 2 „Zl„ ( 19 ) 



The uncertainty equation for u L is determined from 
the minimum and maximum Young's modulus values 
(namely, E min and ii max , respectively) as given below: 



38.330p/ c ^ (L cm -3a L ) 



(22) 



and 



38.330 (p+3<7„)/lzt 



(20) 



and 



38.330p/ c 2 on (L cm +3<7 £ ) 4 



(23) 



where <J p is the estimated one sigma uncertainty of the 
value of p. The uncertainty component, u p , is therefore 
calculated in the same way as before: 



where O l is the estimated one sigma uncertainty of the 
value of L can . Therefore, u L is calculated as follows: 



max mm 



(21) 



E - E 

max min 



(24) 
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The uncertainty equation for u fl . eq is determined from 
the minimum and maximum Young's modulus values 
(namely, E min and E^, respectively) as given below: 



minimum and maximum Young's modulus values 
(namely, E min and is max , respectively) as given below: 



38.330p(/„,„-3a^) 2 C (25) 



38.330p(/ M „-3a /e ) 2 4„ (31) 



and 



and 



38.330p(/ OT +3<T /res ) 2 4, (26 ) 



38.330p(/„,„+3a /e ) 2 4„ ( 32 ) 



where O freq is the one sigma uncertainty of the value of 
f ca „. Therefore, u fi . eq is calculated as follows: 



E - E 

max min 



./''«/ 



(27) 



The uncertainty equation for u p . esol is determined from 
the minimum and maximum Young's modulus values 
(namely, E min and E^, respectively) as given below: 



38.330p|/ ca „-%^ 



(28) 



and 



38.330p 



f f \ 

J r 



f + : 

J can 



L 4 



(29) 



where f resol is the frequency resolution for the given set 
of measurement conditions. Assuming a uniform (i.e., 
rectangular) probability distribution (and assuming 
Up, u l, Ufreq, an d u d equal zero), the value for 



u 



thick J u p J 

E lies between E min and £ M . Therefore, u fresol is calcu- 
lated as follows: 



E. -E_ 



iS 



(30) 



If undamped resonance frequencies (e.g., if the 
measurements were performed in a vacuum) were 
recorded as/ measl ,/ meas2 , and/ meas3 , then u damp is set equal 
to 0.0 Pa. For damped resonance frequencies (i.e., if 
./mea S i,./mea S 2, and /^ were damped measurements), the 
uncertainty equation for u damp is determined from the 



where the calculation for C7 /B is given in Sec. 8.6. 
Therefore, u damp is calculated as follows: 



E - E 

max min 



(33) 



' iitmip 



4.3. Step Height Uncertainty Equations 

In this section, a combined standard uncertainty 
equation is presented for use with step height measure- 
ments (w cSH ). Six sources of uncertainty are identified 
with all other sources of uncertainty considered 
negligible. The six sources of uncertainty are the uncer- 
tainty of the measurement across the length of the step 
(w istep ) where the length is measured perpendicular to 
the edge of the step, the uncertainty of the measurement 
across the width of the step (u Wstep ) where the width is 
measured parallel to the edge of the step, the uncertain- 
ty of the value of the physical step height used for 
calibration (u cert ), the uncertainty of the measurement 
due to the repeatability (w re p eat ), the uncertainty due to 
the amount of drift during the data session (H drift ), and 
the uncertainty of a measurement due to the deviation 
from height linearity of the data scan (wi inear ). As such, 
the combined standard uncertainty equation (as calcu- 
lated in SEMI Test Method MS2 [6] and in Data 
Analysis Sheet SH.l [10]) with six sources of uncer- 
tainty is as follows: 



: V« 



2 2 2 2 

istep U Wstep U cert ~*~ Hepeat 



2 2 

+ *4iff + linear 



(34) 



In determining the combined standard uncertainty, a 
Type B evaluation [17] (i.e., one that uses means other 
than the statistical Type A analysis) is used for each Type 
B source of error, except where noted. Table 4 gives 
sample values for each of these uncertainty components. 
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Table 4. Sample step height uncertainty values 



uncertainty 



Type A or Type B 



value (in |Xm) 



u Lstep 



(using a pkttNXme = 0.0118 urn, 
<V«ivyave = - 0102 li m > 

an d O'roughM' =< 3'roughM' =0 - 00 36lim) 



Type B 



0.011 



u Wstep 



(using stepNxYa = 0.4928 um, 
stepNxft, = 0.4814 um, 
and stepNxYc = 0.4949 urn) 



Type A 



0.0073 



"drift 



(using cert = 9.887 um, 

O cert = 0.083 urn, 

and stepNxY= 0.490 um) 



(using z repeat = 0.024 um 
and Z 6 = 9.879 um) 



(using z drift = 0.016 um 
and Z= 9.887 um) 



(using z pCTC = 1.0%) 



Type B 



Type B 



Type B 



Type B 



0.0041 



0.00034 



0.00023 



0.0028 



"cSH 



3»,,, 



$ 



U, , + Mo, , + U t + U t + U A . f , + li. 

Lstep Wstep cert repeat drift nil 



0.014 



0.041 



This 3w cSH uncertainty is associated with the first repeatability data point plotted in Fig. 7 (for TS1 in quad 1). 



The uncertainty equation for w islep assuming a 
Gaussian distribution is as follows: 



and O piamm . 10 Also, in the above equation, O pUNXwi and 
GptaNYsK are calculated using the following equations: 



' /.step 



[y p/atNXavc ^roughTVX ) """ \y phtNYwc *AoiighiV7 ) ' 

(35) 



" ' platNXn """ ®platNXb """ ® platNX c (36) 



platNX&ve 



and 



where C7 roughMy and C7 roughAT are the surface roughnesses 
of platNX and platNY, respectively, and are calculated 
from the smallest of all the calibrated values obtained 
for fl^Hj and o platNYt , respectively. However, if the 
surfaces of platNX (defined in the List of Symbols), 
platNY, and platNv all have identical compositions, 
then C7 r0llghAft - equals <7 mugbNY , which equals the smallest 
of all the values obtained for O platNXt , (T platNYt , 



" platNY*. OplatNYb "*" ^ platNYc 



platNYsve 



The derivation ofu Lsl is given in Sec. 9.3. 



(37) 



Consult platNrDt in the List of Symbols for the nomenclature 
used in the subscript of (J i atNrDl . 
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The uncertainty equation for u Wstep is determined 
from O mtep , the one sigma standard deviation of the step 
height measurements stepN^, stepN^, and stepN^, 
using the following equation: 



IVsrcp 



Wstep ' 



(38) 



This is a statistical Type A component. 

The uncertainty equation for u cert is determined 
from the minimum and maximum step height values 
(namely, stepN^^ and stepN^^, respectively). The 
uncertainty of the measured step height is assumed to 
scale linearly with height. With 99.7 % confidence, 
assuming a Gaussian distribution (and assuming u Lstep , 
Uwstep, "repeat, "drift, and « linear equal zero), the value for 
\stepNx^ lies between stepNxY min and stepNxy m£i . As 
such, u mrt can be calculated using the following equa- 
tion: 



u 



u 



Lstep , u Wstep , u cert , "drift , 



U, 



U,l 



and w linear equal zero), the value 
for \stepNxA lies between stepNxY min and stepN^^- 
Therefore, M repeat is calculated as follows: 



stepN x 



■ stepN x 



2^ 



which simplifies to the following equation: 



stepNx . 



iSz 



(42) 



(43) 



The uncertainty equation for u^ is determined from 
the minimum and maximum step height values (name- 
ly, stepNxY mm and stepN-^^, respectively) as given 
below: 



cert 



\stepN x 



(39) 



stepN x 



\stepNxY -\stepN x 



2z 



(44) 



where cert is the certified value of the double-sided 
physical step height used for calibration and O cen is the 
one sigma uncertainty of the calibrated physical step 
height. 

The uncertainty equation for M repeat is determined 
from the minimum and maximum step height values 
(namely, stepN^^ and stepN^^, respectively) as 
given below: 



and 



stepNxy^ = IstepN^ -\stepN x 



and 



stepN^^ = IstepN^ + \stepN x 



2z, 



2z, 



(40) 



(41) 



stepN x 



\stepNxY + \stepN x 



2z 



(45) 



where z drift is the calibrated positive difference between 
the averages of the 6 calibration measurements taken 
before and after the data session (at the same location 
on the physical step height used for calibration) and 
where z is the calibrated average of all 12 calibration 
measurements. The uncertainty of the measured step 
height is assumed to scale linearly with height. 
Assuming a uniform distribution (and assuming u Lstep , 
Uwstep, Ucert, "repeat, and u linear equal zero), the value for 
\stepNx^ lies between stepN^^ and stepNxy^. 
Therefore, u dIift is calculated as follows: 



where z repeat is the maximum of two calibrated values; 
one of which is the positive calibrated difference 
between the minimum and maximum values of the six 
calibration measurements taken at a single location on 
the calibration step before the data session and the other 
of which is the positive calibrated difference between 
the minimum and maximum values of the six calibra- 
tion measurements taken at the same location on the 
calibration step after the data session and where z 6 is the 
calibrated average of the six calibration measurements 
from which z repeat is found. The uncertainty of the meas- 
ured step height is assumed to scale linearly with 
height. Assuming a uniform distribution (and assuming 



stepN x 



- stepN x 



Mf ' 2V3 

which simplifies to the following equation: 



drift 



" drift 

iSz 



stepN xy . 



(46) 



(47) 



The uncertainty equation for M lin(;ar is calculated 
from the minimum and maximum step height values 
(namely, stepNxY min and stepNxY mXi , respectively) as 
given: 
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and 



stepN XY mm = \stepN XY \ - \stepN XY \ z perc (48) 



ste P N xy mm = \stepN XY | + \stepN XY \ z perc , ( 49 ) 



where z perc is the maximum percent deviation from 
linearity over the instrument's total scan range, as 
quoted by the instrument manufacturer. The uncertain- 
ty of the measured step height is assumed to scale 
linearly with height. Assuming a uniform distribution, 
u limar can be calculated using the following equation: 



41 



\stepN x 



(50) 



Consult Sec. 9.2 for the uncertainty calculations 
when two step height test structures are used for a step 
height measurement. 



5. Round Robin Results 

The round robin repeatability and reproducibility 
results are given in Sec. 5.1 for Young's modulus 
measurements and in Sec. 5.2 for step height measure- 
ments. 

The repeatability measurements are performed using 
the same test method, in the same laboratory (NIST), 
by the same operator, with the same equipment, in the 
shortest practicable period of time. These measure- 
ments are done on random test structures. 

For the reproducibility measurements, at least six 
independent data sets (each using a different piece of 
equipment or equipment setup) must be obtained 
following the same test method before the results can 
be recorded in the precision and bias statement of a 
SEMI standard test method. These measurements are 
done on random test structures, as described below. 

5.1 Young's Modulus Round Robin Results 

For the Young's modulus portion of the MEMS 
Young's modulus and step height round robin experi- 
ment, both repeatability and reproducibility data were 
taken. 

The repeatability data were taken at one laboratory 
using a dual beam vibrometer (see Sec. 7.1 for specifics 



associated with the vibrometer). Young's modulus 
values were found from 12 different cantilevers 
4 times, with each Young's modulus value determined 
from the average of 3 resonance frequency measure- 
ments. Therefore, 48 Young's modulus values were 
obtained. Of these values, 16 were from 4 different 
cantilevers with L = 200 um, 16 from 4 different 
cantilevers with L = 300 (am, and 16 from 4 different 
cantilevers with L = 400 urn. 

For the reproducibility data, 8 participants were 
identified." Each participant was supplied with a 
Round Robin Test Chip and asked to obtain a Young's 
modulus value from 3 oxide cantilevers with design 
lengths of 200 um, 300 um, and 400 um. (The partici- 
pant could choose to measure any one of five can- 
tilevers of the given length that were available on the 
test chip as long as it passed a visual inspection.) Each 
Young's modulus value was determined from the aver- 
age of three resonance frequency measurements from 
the cantilever as specified in Sec. 3.1, using an instru- 
ment that meets the manufacturer's alignment and 
calibration criteria. Following SEMI MS4 for Young's 
modulus measurements, the measurements were 
recorded on Data Analysis Sheet YM.l [10]. 

The eight participants used a variety of instruments 
(consult Sec. 7.1 for details associated with the instru- 
ments) to obtain Young's modulus. These included a 
single beam vibrometer, a dual beam vibrometer, and a 
stroboscopic interferometer. In addition, thermal exci- 
tation measurements are included for comparison with 
PZT excitation measurements on the same chip. 

The repeatability and reproducibility data for 
Young's modulus and for the combined standard uncer- 
tainty is presented in Table 5 and Table 6, respectively, 
where n indicates the number of calculated Young's 
modulus values. The average of the repeatability or 
reproducibility data (namely E me ) is listed next, 
followed by the 95 % limits for E that is calculated as 
follows: (a) the standard deviations were found, 
(b) these values were multiplied by 2.0 (assuming a 
Gaussian distribution) [17], and (c) the resulting values 
were reported as percents. Below this, the average of 
the combined standard uncertainty [17] values (M cave ) 
and the 95 % limits for u cE are presented. 



The term participant refers to a single data set from a unique com- 
bination of measurement setup and researcher. That is, a single 
researcher equipped with multiple, unique instruments (e.g., a dual 
beam vibrometer and a single beam vibrometer) or different forms of 
excitation (e.g., PZT and thermal excitation) could serve as multiple 
"participants" in this round robin. 
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Table 5. Young's modulus repeatability data (1 participant, 1 laboratory, 1 instrument, 1 chip, 12 different 
cantilevers) 





200 urn 


300 um 


400 nm 


200 


urn to 400 um 




length 


length 


length 




lengths 


n 


16 


16 


16 




48 


£ ave (in GPa) 


59.8 


65.4 


67.5 




64.2 


95 % limits for E 


±1.4% 


± 0.5 % 


± 1.1 % 




± 10.3 % 


"cave ( in G p a) 


2.9 


3.2 


3.3 




3.1 




(4.9 %) 


(4.8 %) 


(4.8 %) 




(4.8 %) 


95 % limits for u cE 


±1.4% 


± 0.5 % 


± 1.1 % 




±10.1 % 



Table 6. Young's modulus reproducibility data (eight participants, five laboratories, seven instruments, four 
chips) 







200 um 


300 um 


400 urn 


200 


um to 400 um 




length 


length 


length 




lengths 


n 


8 


8 


8 




24 


E ave (in GPa) 


58.7 


63.7 


66.0 




62.8 


95 % limits for E 


± 4.4 % 


± 5.5 % 


± 4.4 % 




±11.0% 


"cave ( in GPa) 


2.8 


3.1 


3.2 




3.0 




(4.9 %) 


(4.8 %) 


(4.9 %) 




(4.9 %) 


95 % limits for u cE 


± 4.4 % 


± 5.5 % 


± 4.5 % 




±11.0% 



The Young's modulus round robin results are plotted 
in Fig. 6, where both the repeatability and repro- 
ducibility data are plotted. The average Young's 
modulus value for the repeatability data is specified at 
the top of Fig. 6 along with the average 3u cE uncertain- 
ty bars for this value. 12 These quantities are plotted in 
this figure with both the repeatability and reproducibil- 
ity data. As an observation, all of the reproducibility 
results fall comfortably between the repeatability 
bounds of E me plus or minus 2>u cme . 



Table 3 specifies the value of each of the uncertainty components 
comprising the 3u cE uncertainty bars for the repeatability data point 
corresponding to the first cantilever in Fig. 6 with a length of 
300 um. 



In Fig. 6, the repeatability data are grouped accord- 
ing to the cantilever length with the L = 200 um data 
plotted first, followed by the L = 300 um data, then the 
L - 400 (am data. In like manner with the reproducibil- 
ity data, for each participant, the L = 200 (am data are 
plotted first, followed by the L = 300 um data, then the 
L = 400 um data. The repeatability data and the repro- 
ducibility data both indicate a length dependency. The 
repeatability data in Fig. 6 show a clustering of the data 
at each length. In other words, the 95 % limits for E at 
each length (which are plotted in this figure along with 
E ave for each length) are all less than 1.5 %, which is 
much less than the 10.3 % value (as given in Table 5) 
when all the lengths are considered. This suggests that 
when Young's modulus values extracted by different 
measurement instruments or excitation methods are 
compared, the cantilevers should have the same length. 
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I 
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80 
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100 



Fig. 6. Young's modulus round robin results. 



This length dependency can be due to a number of 
things including debris in the attachment corners of the 
cantilevers to the beam support, which would cause 
larger errors for shorter length cantilevers. This can be 
a topic for future investigation where a) the physical 
form and chemical composition of the cantilever is 
checked to see if it matches the assumptions used in 
the calculations and b) finite element methods are 
used to determine if the length dependency is due 
to the attachment conditions. Therefore, at this point, 
we can only state that, given the existing canti- 
levers, we can only report an "effective" value for 
Young's modulus. 

Round robin participant #1, participant #2, and 
participant #3 took measurements on the same chip 
(chip #1) using a dual beam vibrometer, a single beam 
vibrometer, and a stroboscopic interferometer, respec- 
tively. The results given in Fig. 6 indicate that compa- 
rable results were obtained from these instruments. 

Round robin participant #4, participant #5, and 
participant #6 took data from the same chip (chip #2); 
however, round robin participant #5 used thermal 



excitation to obtain the required data while participant 
#4 and participant #6 used PZT excitation. No signifi- 
cant difference in the results for these measurements is 
seen in Fig. 6. 

No information can be presented on the bias of the 
procedure in the test method for measuring Young's 
modulus because there is not a certified MEMS materi- 
al for this purpose. Many values for Young's modulus 
for various materials have been published with an 
attempt to consolidate this information in [23]. For a 
silicon dioxide film, the Young's modulus values 
reported in [23] range from 46 GPa to 92 GPa. The 
average repeatability value reported in Table 5 of 
64.2 GPa falls comfortably within this range. 

The Young's modulus results are reported as follows 
[17]: Since it can be assumed that the possible estimat- 
ed values are either approximately uniformly distrib- 
uted or Gaussian (as specified in Sec. 4.2) with approx- 
imate standard deviation u cE , the Young's modulus 
value is believed to lie in the interval E ± u cE with a 
level of confidence of approximately 68 % assuming a 
Gaussian distribution. 
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5.2 Step Height Round Robin Results 

For the step height portion of the MEMS Young's 
modulus and step height round robin experiment, both 
repeatability and reproducibility data were taken. 

The repeatability data were taken at one laboratory 
using a stroboscopic interferometer operated in the 
static mode (see Appendix A). Four step height meas- 
urements were taken from each of the four test struc- 
tures in each of the three quads shown in Fig. 1. 
Therefore, 48 step height measurements were obtained 
with 1 step height measurement defined as the average 
of 3 measurements taken from different positions 
somewhat evenly spaced along the step as specified in 
Sec. 3.2. 

For the reproducibility data, seven participants were 
identified. Each participant was supplied with a Round 
Robin Test Chip and asked to obtain the step height 
from any two test structures in the first of the three 
quads of step height test structures. Following SEMI 
MS2 [6] for step height measurements, the raw, uncal- 
ibrated measurements were recorded on Data Analysis 
Sheet SH.l [10]. 

The repeatability and reproducibility data for step 
height and for the combined standard uncertainty is 
presented in Table 7, where n indicates the number of 
step height measurements. The average of the repeat- 
ability or reproducibility data (namely \stepN ABms \) is 
listed next followed by the 95 % limits for \stepN m \ that 
is calculated as follows: (a) the standard deviations 
were found, (b) these values were multiplied by 2.0 
(assuming a Gaussian distribution) [17], and (c) the 
resulting values were reported as percents. Below this, 
the average of the combined standard uncertainty [17] 
values (u cave ) and the 95 % limits for u cSH are 
presented. (It is interesting in comparing the 95 % 



limits for \stepN m \ that the repeatability limits are 
larger than the reproducibility limits. The reason for 
this anomaly is not known.) 

The step height round robin results are plotted in 
Fig. 7 and Fig. 8. In each of these figures, the repeata- 
bility data are plotted first, followed by the results from 
the seven participants. 13 The absolute value of the aver- 
age step height for the repeatability data is specified at 
the top of Fig. 7 and Fig. 8 along with the 3w cave uncer- 
tainty bars for this value, as obtained or derived from 
Table 7. 14 These quantities are plotted in each figure 
with both the repeatability and reproducibility data. As 
an observation, all of the reproducibility results fall 
comfortably between {stepN^^l pl us or minus 3u cave as 
obtained from the repeatability results. 

Figure 7 groups the repeatability results by quad 
number with the results from quad "1" plotted first, fol- 
lowed by the results from quad "2," then the results 
from quad "3." The results within each quad are 
grouped according to test structure number 15 with the 
results from test structure "1" plotted first, followed by 
the results from test structure "2," etc. The average step 
height value and the 95 % limits for this value for 
each quad are given at the bottom of Fig. 7 and also in 
Table 8. These results reveal comparable values for the 
step height measurements and comparable values for 
the 95 % limits. This implies there are no discernable 
variations in the step height value between neighboring 
quads. 

Figure 8 groups the repeatability results by test struc- 
ture number with the results from test structure "1" 
(TS1) plotted first, followed by the results from test 
structure "2" (TS2), followed by the results from test 
structure "3" (TS3), then test structure "4" (TS4). The 
results for each test structure number are grouped 
according to quad with the results from quad "1" 



Table 7. Step height measurement results 





Repeatability 


Reproducibility 




results 


results 


n 


48 


14 


\stepN mm<: \ (in um) 


0.477 


0.481 


95 % limits for \stepN AB \ 


+ 7.9 % 


+ 6.2 % 


"cave ( in um ) 


0.014 


0.014 




(3.0 %) 


(3.0 %) 


95 % limits for w cSH 


+ 70 % 


+ 80 % 



Participant #2 provided stylus profilometer results, and that 
instrument did not provide standard deviation values; however it did 
provide average roughness values. Although this test method uses 
standard deviation values, in this analysis for that laboratory average 
roughness values were inserted into the data sheet for analysis as 
opposed to standard deviation values. 

Table 4 specifies the value of each of the uncertainty components 
comprising the 3w cSH uncertainty bars for the first repeatability data 
point plotted in Fig. 7 (for TS1 in quad 1). 

The upper left hand step height test structure in a quad, such as 
shown in Fig. 4, is called test structure "1" and it has a 0° orientation. 
Test structure "2" (the upper right test structure) has a 270° orienta- 
tion. Test structure "3" (the bottom right test structure) has a 
180° orientation. And test structure "4" (the bottom left test structure) 
has a 90° orientation. 
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Fig. 7. Step height round robin results with the repeatability results grouped according to quad 
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Fig. 8. Step height round robin results with the repeatability results grouped according to test structure number. 
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plotted first, followed by the results from quad "2," 
then the results from quad "3." As in Fig. 7, the average 
step height value and the 95 % limits for this value for 
each test structure are given at the bottom of Fig. 8 and 
also in Table 9. These results reveal that TS1 and TS3 
(which are rotated ± 90° with respect to TS2 and TS4) 
have comparable 95 % limits as do TS2 and TS4; 
however the 95 % limits for TS1 and TS3 are slightly 
less than the 95 % limits for TS2 and TS4 when they 
should be comparable. There are also more variations 
in the average step height value between rotated test 
structures (as shown in Fig. 8 and Table 9) than varia- 
tions in this value between quads (as shown in Fig. 7 
and Table 8). 

Table 8. Step height repeatability data grouped by quad 







Ql 


Q2 


Q3 


11 


16 


16 


16 


l-SteP^ABavel 


0.479 


0.473 


0.478 


(in urn) 








95 % limits for 


+ 8.4 % 


± 7.2 % 


+ 8.3 % 


\stepN m \ 








"cave ( in 1™) 


0.015 


0.015 


0.013 




(3.1 %) 


(3.2 %) 


(2.8 %) 


95 % limits for 


+ 200 % 


± 200 % 


+ 200 % 


"cSH 









platform) and the data jagged. Therefore, the precise 
selection of the analysis regions (including the number 
of data points within these regions) affects the standard 
deviations obtained. An averaging capability incorpo- 
rated in most non-contact instruments can have the 
effect of smoothing the data; however, a more compre- 
hensive determination of the length and width varia- 
tions may be necessary when dealing with tilt. 
Repeatability might also be improved by calculating 
the step height from fitted straight lines as described for 
NIST step height calibrations [24] and outlined in 
ASTM E2530 [25]. As given in [24], "For step height 
measurements, one of several algorithms may be used. 
For single-sided steps, a straight line is fitted by the 
method of least squares to each side of the step transi- 
tion, and the height is calculated from the relative posi- 
tion of these two lines extrapolated to the step edge." 

The interferometer or comparable instrument is cali- 
brated in the out-of-plane z-direction. If the calibration 
is not done, a bias to the measurements is expected. The 
direction and degree of the resulting bias is different for 
each magnification of each instrument. As such, cali- 
bration of the interferometer or comparable instrument 
is considered mandatory for step height measurements. 

The step height results are reported as follows [17]: 
Since it can be assumed that the possible estimated 
values are either approximately uniformly distributed 
or Gaussian (as specified in Sec. 4.3) with approximate 
standard deviation u cSH , the step height is believed 
to lie in the interval stepNxy+ u cSH with a level of con- 
fidence of approximately 68 % assuming a Gaussian 
distribution. 



Table 9. Step height repeatability data grouped by test structure 



TS1 



TS2 



TS3 



TS4 



n 


12 


12 


12 


12 


\stepN mw ^ 


0.486 


0.469 


0.474 


0.478 


(in nm) 










95 % limits 


+ 6.4 % 


± 8.7 % 


± 7.0 % 


+ 8.5 % 


for \stepNp^\ 










"cave ( in nm) 


0.014 


0.015 


0.011 


0.018 




(2.8 %) 


(3.1 %) 


(2.4 %) 


(3.7 %) 



95% limits +28.3% +71.6% ±27.7% ±81.2% 



for u 



cSH 



The platform surfaces involved in the step were not 
ideal surfaces. Oftentimes they were tilted (even though 
the data were leveled with respect to the reference 



6. Conclusions 

The technical basis for the two SEMI standard test 
methods on Young's modulus measurements [5] and 
step height measurements [6] was presented, along 
with the data obtained from the MEMS Young's modu- 
lus and step height round robin experiment. These data 
were incorporated into the precision and bias state- 
ments for the applicable SEMI standards. 

For Young's modulus measurements, data obtained 
from a single beam vibrometer, a dual beam vibrome- 
ter, and a stroboscopic interferometer yielded compara- 
ble results. Also, PZT excitation and thermal excitation 
yielded comparable results. The repeatability data and 
the reproducibility data both indicate a length depend- 
ency. In other words, the 95 % limits for E at each 
length are all less than ±1.5 %, which is much less than 
the + 10.3 % value when all the lengths from 200 urn to 
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400 um are considered. This could be due to debris in 
the attachment corners of the cantilevers to the beam 
support that would cause larger errors for shorter length 
cantilevers. Additional research can ascertain if the 
physical form and chemical composition of the 
cantilever matches the assumptions used in the calcula- 
tions and if the length dependency is due to the attach- 
ment conditions. Therefore, at this point, we can only 
state that, given the existing cantilevers, we can only 
report an "effective" value for Young's modulus. Also, 
when Young's modulus values extracted by different 
measurement instruments or excitation methods are 
compared, the cantilevers should have the same length. 

For step height measurements, it is currently not 
understood why the repeatability limits (± 7.9 %) are 
larger than the reproducibility limits (± 6.2 %). The 
platform surfaces involved in the step are not ideal 
surfaces, such that the precise selection of the analysis 
regions (including the number of data points within 
these regions) affects the standard deviations obtained. 
Also, the repeatability data indicate there are no 
discernable variations in the step height value between 
neighboring quads. And, there are more variations in 
the average step height value between rotated test 
structures than variations in this value between quads. 
For one-sided steps, the large uncertainties can be 
improved by calculating the step height from fitted 
straight lines (using the method of least squares) with 
the height calculated from the relative position of the 
two lines extrapolated to the step edge. This approach 
is typically used for step height measurements and 
calibrations [24,25]. 

It is expected that these standards will be instrumen- 
tal in reducing the interlaboratory differences in the 
parametric measurements. The following guidelines 
should be of value to the MEMS industry, when 
communicating data obtained using these SEMI 
standard test methods: 

(a) To record the data and perform the calculations, 
use the data analysis sheets which are accessible 
via the NIST SED MEMS Calculator website 
[10], 

(b) Make sure the data have passed the verification 
check given at the bottom of the data analysis 
sheets, and 

(c) Ask if all the verification checks were passed, 
when communicating with others concerning 
SEMI standard data. 



Consult Ref [10] for information associated with the 
MNT 5-in-l SRM that enables companies to compare 
their in-house measurements taken on this SRM with 
NIST measurements thereby validating their use of the 
documentary standards. 



7. Appendix A. 

Instrument Specifications 

An optical vibrometer, stroboscopic interferometer, 
or comparable instrument is required for Young's 
modulus measurements, as specified in Sec. 7.1. For 
step height measurements, an optical interferometer or 
comparable instrument is used as specified in Sec. 7.2. 

7.1 Instrument Specifications for Young's 
Modulus Measurements 

For Young's modulus measurements, an optical 
vibrometer, stroboscopic interferometer, or an instru- 
ment comparable to one of these is required that is 
capable of non-contact measurements of surface 
motion. This section briefly describes the operation and 
specifications for a typical single beam laser vibro- 
meter, a dual beam laser vibrometer, and a stroboscop- 
ic interferometer. The specifications can be applied to 
comparable instruments, if appropriate. 

For a single beam laser vibrometer, a typical 
schematic is given in Fig. 9. A signal generator pro- 
vides an excitation signal, which excites the sample via 
a piezoelectric transducer (PZT). The measurement 
beam is positioned to a scan point on the sample (by 
means of mirrors) and is reflected back. The reflected 
laser light interferes with the reference beam at the 
beam splitter (BS). A photodetector (PD) records the 
interference signal, converting it into an electrical 
signal. The frequency difference between the beams is 
proportional to the instantaneous velocity of the vibra- 
tion parallel to the measurement beam. (The Bragg cell 
is instrumental in determining the sign of the velocity.) 
The velocity decoder provides a voltage proportional to 
the instantaneous velocity. 

A dual beam laser vibrometer incorporates two 
beams. The measurement beam is positioned to a scan 
point on the sample (for example, positioned near the 
tip of a cantilever). The reference beam emanates from 
the beam splitter (BS) shown in Fig. 9 and is positioned 
to a second point on the sample (for example, posi- 
tioned on the support region at the base of the can- 
tilever). The two scattered beams optically combine at 
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Fig. 9. Schematic of a typical setup for a single beam laser vibrometer. (PBS indicates a polarizing beam 
splitter, BS indicates a beam splitter, P indicates a prism, and PD indicates a photodetector.) 



the beam splitter where the reference beam is used to 
directly eliminate any movement of the sample also 
experienced by the measurement beam. 

For a stroboscopioc interferometer, a simplified 
schematic of a typical setup is shown in Fig. 10(a). 
When operated in the static mode, the interferometer is 
used to determine surface profiles. The incident light 
travels through the microscope objective to the beam 
splitter. Half of the light travels to the sample surface 
and then back to the beam splitter. The other half is 
reflected to a reference surface and then back to the 
beam splitter. These two paths of light recombine at the 
beam splitter to form interference light fringes. As the 
interferometer scans downward, an intensity envelope 



incorporating these fringes is determined by the soft- 
ware [see Fig. 10(b)]. The peak contrast of the fringes, 
phase, or both are used in determining the sample 
height at that pixel location. The surface profile is 
found by collecting sample height data for each pixel 
within the field of view (FOV). When operated in the 
dynamic mode, the incident light is strobed at the same 
frequency as that used to actuate the device. The 
sample is actuated, for example, after securing it to the 
top of a PZT. The phase, frequency, and drive signal to 
the strobe and PZT are varied, performing a downward 
scan as is done for static devices at each combination to 
obtain successive 3D images as the sample cycles 
through its range of motion. 
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Fig. 10. For a typical stroboscopic interferometer (a) a schematic and (b) an intensity envelope used to obtain 
a pixel's sample height. 
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Specifications for the (previous) above instruments or a 
comparable instrument are as follows: 

1. The microscope objective or objectives should 
be chosen to allow for sufficient resolution of the 
cantilever or fixed-fixed beam and a portion of the 
surrounding sample. The objective(s) should have a 
FOV that can encompass at least half of the length of 
the cantilever or fixed-fixed beam being measured. 
Typically, a 4x and a 20x objective will suffice. 

2. The signal generator should be able to produce a 
waveform function (such as a periodic chirp function 16 
or a sine wave function 17 ) if applicable, such that from 
its use, a reproducible resonance frequency can be 
obtained and good 3-D oscillating images can be 
obtained such that it is obvious by inspection that the 
beam is in resonance. 

3. The instrument shall be capable of producing a 
magnitude versus frequency plot from which the 
resonance frequency can be obtained. 

4. The instrument should be capable of obtaining 3-D 
images of beam oscillations in order to ascertain if the 
correct frequency peak has been chosen as the beam's 
resonance frequency. 

5. An estimate for the maximum frequency of the 
instrument needed for a resonating cantilever, f cantni „ is 
at least the value calculated using the following 
equation: 18 



J ca 



e.,j z 



|38.330pr 



(51) 



The periodic chirp function enables quick results without averag- 
ing. For the periodic chirp function, sinusoidal signals (within the 
selected frequency range and of the same amplitude) are emitted 
simultaneously for all fast Fourier transform (FFT) lines. The period- 
ic chirp function is periodic within the time window, and the phases 
are adapted to maximize the energy of the resulting signal. 

The periodic chirp function produces a reproducible resonance 
frequency. A sine wave sweep function produces a resonance 
frequency that can be affected by the direction of the sweep if there 
is not a sufficient amount of time between measurements. 

By inserting the inputs into the correct locations on the appropri- 
ate NIST Web page [10], the given calculation can be performed 
on-line in a matter of seconds. 



as derived in Sec. 8.1. An estimate for the maximum 
frequency of the instrument needed for a resonating 
fixed-fixed beam,f f/bi „ ithi , is at least the value calculated 
using the following equation: 



./; 



£„„r 



lihiiiiiiii 



as derived in Sec. 8.5. 



1 0.946pZ 4 



(52) 



///. 



6. An instrument that can make differential measure- 
ments (e.g., with the use of two laser beams) is recom- 
mended for use with fixed-fixed beams. It is also 
recommended for use with cantilevers, especially for 
estimated resonance frequencies less than 10 kHz and 
also if the value for p djff as calculated in the following 
equation is greater than or equal to 2 %: 



Pdiff 



l-.l 



1 



\ 



4Q 2 



100%. 



(53) 



For a cantilever, the Q-factor, Q, in the above equation 
can be estimated using the following equation: 



Q- 



w m 4K 

2Afi 



( t 




L 



(54) 



7.2 



Instrument Specifications for Step Height 
Measurements 



For step height measurements, an optical interfero- 
metric microscope or comparable instrument is used 
that is capable of obtaining topographical 2-D data 
traces. (The stroboscopic interferometer operated in the 
static mode, as described in Sec. 7.1, can be used for 
these measurements.) Figure 11 is a schematic of a typ- 
ical optical interferometric microscope that uses the 
method of coherence scanning interferometry [26], also 
called vertical scanning interferometry or scanning 
white light interferometry, for these measurements. 
However, any calibrated topography measuring instru- 
ment that has pixel-to-pixel spacings or sampling 
intervals as specified in Table 10 and that is 
capable of performing the test procedure with a vertical 
resolution less than 1 nm is permitted. The optical 
interferometric microscope or comparable instrument 
must be capable of measuring step heights to at least 5 
urn higher than the step heights to be measured and 
must be capable of extracting standard deviations. 
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Fig. 11. Schematic of an optical interferometric microscope operating in the Mirau configuration where the 
beam splitter and the reference surface are between the microscope objective and the sample. 



Table 10. Interferometer pixel-to-pixel spacing requirements 



Magnification, X 



Pixel-to-pixel spacing, um 



5 

10 
20 
40 



< 1.57 
<0.83 
<0.39 

< 0.21 

< 0.11 



This table does not include magnifications at or less than 2.5 X for 
optical interferometry because the pixel-to-pixel spacings will be too 
large for this work and the possible introduction of a second set of 
interferometric fringes in the data set at these magnifications can 
adversely affect the data. Therefore, magnifications at or less than 
2.5 x shall not be used with optical interferometry. 



Appendix B. 

Supplemental Material for Young's 

Modulus Measurements 



deriviation of the Young's modulus equation from a 
resonating fixed-fixed beam assuming simply support- 
ed boundary conditions is given in Sec. 8.4, and the 
derivation of the Young's modulus equation from a 
resonating fixed-fixed beam assuming clamped- 
clamped boundary conditions is given in Sec. 8.5. 
Lastly, Sec. 8.6 presents the calculation of Cyg as 
discussed in Sec. 4.2. 

8.1 Derivation of Young's Modulus Equation 
From a Resonating Cantilever 

In this section, the Young's modulus equation is pre- 
sented as derived, in a manner similar to that presented 
in [12], from measurements taken on a resonating 
cantilever. Clamped-free boundary conditions are 
assumed, as given by the following equations: 



(*L =0: 



(55) 



This appendix contains supplemental material for 
Young's modulus measurements for reference only. In 
Sec. 8.1, the derivation of the Young's modulus equa- 
tion from a resonating cantilever is presented. In 
Sec. 8.2, the residual stress and stress gradient equa- 
tions are given. The equations for Young's modulus 
measurements as obtained from a fixed-fixed beam are 
given in Sec. 8.3 followed by two derivations. The 



and 



(dX^ 



. dx . 



U 2 x^ 



dx 2 
V J 



:0, 



(56) 



(57) 
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<d'x^ 



dx' 



= 0. 



(58) 



where the cantilever attachment point is at x = um, 
where the tip of the cantilever is at x = L, where L is the 
length of the cantilever, and where X is the normal 
function [12]. 

The equation of motion for all free transverse vibra- 
tions is given by the following equation [12]: 



( d'X^ 



dx 



f n A\ 



pA 
K EI J 



0) 2 X = . 



(59) 



where E, I, p, A, and (0 are the Young's modulus, the 
area moment of inertia about the neutral axis, the 
density, the area, and the angular frequency, respective- 
ly. To help solve this fourth-order differential equation, 
the following notation is introduced: 



0) = ak 



(60) 



such that the equation of motion can be given by the 
following equation: 



( d'X^ 



dx 



-k 4 X = 0. 



(61) 



A solution to the above equation yields: 



X = C X (cos kx + cosh kx) + C 2 (cos kx -cosh kx) 
+ C 3 (sin kx + sinh kx) + C 4 (sin kx -sinh kx) . 



(62) 



[This can be demonstrated by differentiating Eq. (62) 
four times and plugging the result into Eq. (61) along 
with Eq. (62) to achieve an identity.] Applying the four 
boundary conditions yields the following equation: 



(cos 2 kL + sin 2 kL\+ (cosh 2 /fcZ-sinh 2 kL} 



-2 cos kL cosh kL . 



(63) 



However, since 



and 



cos x + sin x = l 



cosh 2 x-sinh 2 x = l , 



(64) 



(65) 



Equation (63) can be rewritten as 
cos k,L cosh k,L = -1 



(66) 



where k.L ~ 






n with z =1,2,3...°° 



or 



k t L= 1.875, 4.694, 7.855,. 



The frequency equation is given by: 
ft) ak. 



f, 



In In 



where 



Therefore, for k x L = 1.875, we have 



a k?L 2 a f 1.875^ 
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or rewritten 
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where t is the thickness so 
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and 
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38.330pZ 4 



(67) 



(68) 



(69) 
(70) 

(71) 
(72) 



(73) 



(74) 



stated already as Eq. (1) in Sec. 3.1 and Eq. (51) in 
Sec. 7.1. Soving for E produces the following equation: 



38.330p/, 2 Z 4 



(75) 



or Eq. (7) in Sec. 3.1. 
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8.2 Residual Stress and Stress Gradient Equations 

Given the Young's modulus value, E, obtained from 
Eq. (7) in Sec. 3.1, the residual stress, C7 r , and the stress 
gradient, o g , can be found from Eqs. (76) and (80) 
presented in this section. 

To calculate residual stress, the residual strain of the 
thin film layer, e r , must also be known. This value and 
its combined standard uncertainty value, u cer , are found 
from measurements of a fixed- fixed beam test structure 
comprised of that layer using ASTM E2245. (ASTM 
E2245 can be consulted to determine if the fixed-fixed 
beam used for this residual strain measurement can also 
be used to find /a, and E in Sec. 8.3. However, it is 
recommended that the Young's modulus value be 
obtained from the average resonance frequency of a 
cantilever as specified in Sec. 3.1 due to a lower result- 
ing value for u cE .) Then, the residual stress is calculat- 
ed using the following equation: 



V(o,) 



(79) 



where u cer was found in ASTM E2245 for residual 
strain. 

To calculate stress gradient, the strain gradient of the 
thin film layer, s g , must also be known. This value and 
its combined standard uncertainty value, u csg , are found 
from measurements of a cantilever test structure 
comprised of that layer using ASTM E2246. (ASTM 
E2246 can be consulted to determine if the cantilever 
used for this strain gradient measurement can also be 
used to find f can and E in Sec. 3.1.) Then, the stress 
gradient, o g , is calculated using the following equation: 



■Es, 



(80) 



The combined standard uncertainty for stress gradient, 
w CCTg , is then calculated using the following equation: 



Ee, 



(76) 



The combined standard uncertainty for residual 
stress, u car , is then calculated using the folowing 
equation: 



v^5 



. , +u , , 

Br) er(ar) 



(77) 



where u E , ar) and u er(ar) are defined in the List of 
Symbols. In determining the combined standard uncer- 
tainty, a Type B evaluation [17] (i.e., one that uses 
means other than the statistical Type A analysis) is used 
for each source of error. 

The uncertainty equation for u E(ar} in the above equa- 
tion is determined from the minimum and maximum 
residual stress values (namely, a min and (7,™*, respec- 
tively). Assuming a Gaussian distribution, u E(ar) is 
calculated using the following equation: 



*E(ar) ■ 



(78) 



where u cE is the combined standard uncertainty of the 
Young's modulus measurement as given in Sec. 4.2 
(or u E from Sec. 8.3 can be used although this is not 
recommended). 

The uncertainty equation for u er(ar) in Eq. (77) is 
determined from the minimum and maximum residual 
stress values (namely, O^ and ff^, respectively). 
Assuming a Gaussian distribution, u ei(C7r) is calculated 
using the following equation: 



■■ft 



2 2 

E("g) sg(og) 



(81) 



where u E((Jg} and u sg(ag) are defined in the List of 
Symbols. In determining the combined standard uncer- 
tainty, a Type B evaluation [ 1 7] (i.e., one that uses 
means other than the statistical Type A analysis) is used 
for each source of error. 

The uncertainty equation for u E{ag) in the above equa- 
tion is determined from the minimum and maximum 
stress gradient values (namely, a min and C7 gmax , respec- 
tively). Assuming a Gaussian distribution, u E(ag) is 
calculated using the following equation: 



**(<*) ■ 



(82) 



where u cE is the combined standard uncertainty of the 
Young's modulus measurement as given in Sec. 4.2 (or 
u E from Sec. 8.3 can be used although this is not recom- 
mended). 

The uncertainty equation for u sg(ag) in Eq. (81) is 
determined from the minimum and maximum stress 
gradient values (namely, (7 gmin and O^, respectively). 



Assuming a Gaussian distribution, u sg{ag) 
using the following equation: 



*«(CTg) 



: U c, S E : 



is calculated 



(83) 



where u csg was found in ASTM E2246 for strain 



gradient. 
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8.3 Young's Modulus Measurements From 
Fixed-Fixed Beams 

The Young's modulus of a single layer can be 
obtained from resonance frequency measurements 
from a fixed-fixed beam comprised of the same materi- 
al (although the use of a cantilever is recommended as 
described in Sec. 3.1 due to a lower resulting combined 
standard uncertainty value). 

To determine an estimate for the resonance frequen- 
cy of a fixed-fixed beam, the lower and upper bounds 
for the resonance frequency (namely, fgumtio and 
fffl, in uu, respectively) are calculated using the following 
equations (as derived in Sec. 8.4 and Sec. 8.5, respec- 
tively): 19 



.///, 



E. ./ 

mil 



|4.864pZ; 



(84) 



Given these two values, E simple and E clmnped , the average 
Young's modulus is calculated using the following 
equation: 



simple 



+ E 



'clamped 



(88) 



and recorded as the Young's modulus value, E, for the 
thin film layer. 

The standard uncertainty, u E , for the Young's modu- 
lus of the film is determined from a Type B evaluation 
[17] using the calculated minimum and maximum 
Young's modulus values (namely, E dampei and E simple , 
respectively). With 99.7 % confidence, assuming a 
Gaussian distribution, the value of E lies between 
E damped and ^simple ■ Therefore, u E is calculated as follows: 



and 



simple 



^clamped 



3.918p^4 (89) 
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(85) 
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1 0.946pZ 4 



where E Mt is the initial estimate for the Young's modu- 
lus value of the thin film layer, t is the thickness of the 
suspended layer, p is the density, and L l]b is the 
suspsended length of the fixed-fixed beam. 

Measurements are taken which encompass the aver- 
age offjpMto waAfjgbMM and an excitation-magnitude ver- 
sus frequency plot is obtained from which the reso- 
nance frequency is found. For a given fixed-fixed 
beam, given the average of at least three resonance 
frequency measurements,^,, two values of Young's 
modulus are calculated; one assuming simply support- 
ed boundary conditions at both supports, E simpk , (as 
derived in Sec. 8.4) and given below: 



48/r 2 pf m r m 4.864py* £ 



.//'' n.lh 



J simple 



(86) 



8.4 Derivation of Young's Modulus Equation 
From a Resonating Fixed-Fixed Beam 
Assuming Simply Supported Boundary 
Conditions 

In this section, the Young's modulus equation is 
derived, in a manner similar to that presented in [12], 
from resonance frequency measurements taken on a 
fixed-fixed beam. Simply supported boundary condi- 
tions are assumed for both supports, as given by the 
following equations: 



(*Lo=°> 



U 2 X^ 



dx 



(90) 



(91) 



and one assuming clamped-clamped boundary condi- 
tions, E clamped , (as derived in Sec. 8.5) and given below: 



and 
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(92) 



48/r 2 pflA, 0.946 pflC 
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4.730 4 



(87) 



(d 2 X^ 



dx' 



= 0, 



(93) 



By inserting the inputs into the correct locations on the appropri- 
ate NIST SED Web page (http://www.eeel.nist.gov/812/test-struc- 
tures/MEMSCalculator.htm), the calculation can be performed on- 
line in a matter of seconds. 



where the fixed-fixed beam attachment points are at 
x = (am and at x = L, where L is the length of the fixed- 
fixed beam, and where X is the normal function [12]. 
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The equation of motion for all free transverse vibra- 
tions is given by the following equation [12]: 



Therefore, with k,L = ffwe have 



U'X^ 



dx 
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where E,I,p,A, and CO are the Young's modulus, the 
area moment of inertia about the neutral axis, the 
density, the area, and the angular frequency, respective- 
ly. To help solve this fourth-order differential equation, 
the following notation is introduced: 



K EI 
h ~ 2L 1 \pA 



(102) 



ft)= ak 



(95) 



But 



such that the equation of motion can be given by the 
following equation: 



I = — At 1 

12 
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dx 4 
v J 



■rz=o. 
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A solution to the above equation yields: 



or rewritten 



II t 



'A Vl2' 



(104) 



X = C, (cos kx + cosh kx) + C 2 (cos kx -cosh kx) /ny\ 
+ C 3 (sin kx + sinh kx) + C 4 (sin kx -sinh kx) . 



where t is the thickness so 



[This can be demonstrated by differentiating Eq. (97) 
four times and plugging the result into Eq. (96) along 
with Eq. (97) to achieve an identity.] Applying the four 
boundary conditions yields the following equation: 



/, 



n t E n Et 
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sin k t L = , 

IK 

where k. = — with i = 1,2,3,..., °°. 



The frequency equation is given by: 



(98) 



and 
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(99) 



[see Eq. (84) in Sec. 8.3]. Solving for E produces the 
following: 



where 



4.864P/, 2 ! 4 



(107) 



(100) 



[see Eq. (86) in Sec. 8.3]. 



333 



Volume 115, Number 5, September-October 2010 

Journal of Research of the National Institute of Standards and Technology 



8.5. Derivation of Young's Modulus Equation 
From a Resonating Fixed-Fixed Beam 
Assuming Clamped-Clamped Boundary 
Conditions 

In this section, the Young's modulus equation is 
derived, in a manner similar to that presented in [12], 
from resonance frequency measurements taken on a 
fixed-fixed beam. Clamped-clamped boundary condi- 
tions are assumed, as given by the following equations: 



(*Lo 



dX^ 

dx 



:0, 



(108) 



(109) 



A solution to the above equation yields: 



X = C, (cos kx + cosh kx) + C 2 (cos Ax -cosh kx) 
+ C 3 (sin kx + sinh kx) + C 4 (sin kx -sinh kx) . 



(115) 



[This can be demonstrated by differentiating Eq. (115) 
four times and plugging the result into Eq. (114) along 
with Eq. (115) to achieve an identity.] Applying 
the four boundary conditions yields the following 
equation: 



(cos 2 kL + sin 2 kL )+ (cosh 2 kL - sinh 2 klA 



2 cos kL cosh kL. 



However, since 



(116) 



and 
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(110) 



(111) 



where the fixed-fixed beam attachment points are at 
x = urn and at x = L, where L is the length of the fixed- 
fixed beam, and where X is the normal function [12]. 

The equation of motion for all free transverse vibra- 
tions is given by the following equation [12]: 



( d"X^ 



dx 
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(112) 



where E, I, p, A, and ft) are the Young's modulus, the 
area moment of inertia about the neutral axis, the 
density, the area, and the angular frequency, respective- 
ly. To help solve this fourth-order differential equation, 
the following notation is introduced: 



ft) = ak 



(113) 



cos 2 x + sin x = 1 



and 



cosh 2 x-sinh 2 x = l 
Eq. (116) can be rewritten as 



cos k L cosh k L =1 , 



/ i \ 



where k,L ~ 



i + - 



or k t L = 4.730, 7.853, 10.996... . 

The frequency equation is given by: 
ft), ak, 
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where 



Therefore, for k t L = 4.730, we have 
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k with i =1,2,3,..., °° 
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such that the equation of motion can be given by the 
following equation: 



< d 4 X^ 



dx 



-k 4 X = 0. 
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or 
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where O w is the estimated one sigma uncertainty of the 
value of W mn and f dampedave is the average of the three 
damped resonance frequency values (f damped i,f damped2 , and 
/damped})- With 99.7 % confidence, assuming a Gaussian 
distribution (and assuming u L(f) , u (f) , u 



(f)> u p(f)> u tff), "E(f) 
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fundampedmm anc * fundampedmax- Therefore, u W (pis calculated as 
follows: 



''/' If I 



equal zero), the value for f u 



undamped 



f- I Et 



(127) 



/;, 



undamped r 



Ju, 



undamped i 



*W(f) 



(132) 



which agrees with Eq. (52) in Sec. 7.1 and Eq. (85) in 
Sec. 8.3. Solving for E produces the following: 



0.946p/ 1 2 i 4 



which agrees with Eq. (87) in Sec. 8.3. 



(128) 



8.6 Calculation of <x 



7!? 



To obtain a fQ (as given in Sec. 4.2), the following 
equation is used: 



— r~ 2 ~i_ ~ 2 ~z ~ 2 ~z ~ 2 ~z 2 a 2 

°fQ ~ \ U W(f) + U L(f) + U p(f) + U l(f) + U E(f) + U p(f) » 



(129) 



where u w(f) , u L(f) , u p(f) , u t(f) , u E(f) , and u fl(f) are found 
below. 

The uncertainty equation for u w(f) is determined from 
a calculation of the potential variation in the undamped 
resonance frequencies due to the uncertainty in the 
cantilever width, W can , as given below: 



/;, 



./; 



tla/ufiL'Ju; 



undamped min 



The uncertainty equation for u L(f) is determined from 
the minimum and maximum undamped resonance 
frequencies (namely, f mdampedmm and f mdampedtwil , respec- 
tively) as given below: 



./„ 



/„ 



dampedave 



undamped min 



KanJzp 



2Aji 



(A»-3ffJ 



(133) 



./;, 



./;, 



dampedave 



ndamped max 



1 



W cm ^Ep{ t 



2Ap. 



(L cm + 3a L ) 



(134) 



where O l is the estimated one sigma uncertainty of the 
value of L ca „. With 99.7% confidence, assuming a 
Gaussian distribution (and assuming u w(f) , u p(f) , u t(f) , 
u E(f) , and u p(f) equal zero), the value for f mJamped lies 
between f u 
calculated as follows 



undampedmin an " Jundampedmxi- IbeiCIOTe, U^y) IS 



(Km+2°wHEP 

2Ajl 



v 4™ j 



(130) 



.' it 



undamped max Jundamped 



Jum 



*£(/) 



(135) 
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The uncertainty equation for u p(f) is determined from 
the minimum and maximum undamped resonance 
frequencies (namely, f undampedmm and f undampedm& , respec- 
tively) as given below: 



/„ 



fa 



dampedave 



undamped min 



1 



W^E{ P + 7,a p ) 
2Ap 



f j_ 

L 

V cm ) 



(136) 



where a thick is the one sigma uncertainty of the value of 
t, which is found using the electro-physical technique 
[18]. Data Sheet T.l [10] can be used to calculate O thick . 
With 99.7 % confidence, assuming a Gaussian distribu- 
tion (and assuming u w(f) , u L(f) , u p(f) , u E(f) , and u p(f) 
equal zero), the value for f undamped lies between 
fundampedmin an d fundampednm. ■ Therefore, u uf) is calculated as 
follows: 



J u 



undamped i 



■ /' 

■ nil 



undamped i 



-«/) 



(141) 



/; 



/;, 



tlmnpci/a) 



undamped max 



1 






< t * 



(137) 



where O p is the estimated one sigma uncertainty of the 

value of p. With 99.7 % confidence, assuming a 

Gaussian distribution (and assuming u 

u E(f) , and u p(f) equal zero), the value for /„„,„„„.., 

between f u 

calculated as follows: 



W(f)> u L(f)i w t(f)> 

lies 

imdampedmm an U J undampednux ■ * fieretOre, Upff) IS 



/;, 



undamped I 



-' u. 



undamped min 



(138) 



p(f) 



The uncertainty equation for u t(f) is determined from 
the minimum and maximum undamped resonance 
frequencies (namely, f undampedmm and f undampedmm , respec- 
tively) as given below: 

/_ J dampedave 

undamped min 



./;, 




undamped max 



W^yfEpftSa. ' 



24,u 



(140) 



The uncertainty equation for u E(f) is determined from 
the minimum and maximum undamped resonance 
frequencies (namely, f mdampedmm and f undampedmm , respec- 
tively) as given below: 



/; 



/; 



dampedave 



undamped min 



1 



W m „4(E+l>u cE0 ) P 



< t * 



2Ap. 



L 



142) 



/;, 



/;, 



dampedave 



undamped max 



1 



W m n^(E-3u cE0 )p 



2<\H 



L 



(143) 



where u cE0 is given by: 



\~2 , ~. ~. 2 , 2 

U r „ = Ju, u . , + U„ + U, + U, + U, , 

c£0 V Thick p L Jreq jresol 



(144) 



With 99.7 % confidence, assuming a Gaussian distribu- 



(139) tion (and assuming u 



. u p(f) ,u t(f) , and u p(f) equal 



W(f) i u L(f) i 

zero), the value for f undmnped lies between 
fundampedmm an d fundampedmux ■ Therefore, u E( f) is calculated as 
follows: 



./„ 



undamped r 



■ 'in 



undamped i 



'£(/) 



(145) 



336 



Volume 115, Number 5, September-October 2010 

Journal of Research of the National Institute of Standards and Technology 



The uncertainty equation for u ll(f) is determined from 
the minimum and maximum undamped resonance 
frequencies (namely, f undampedmm and f umiampedmm , respec- 
tively) as given below: 

/__ J dampedave 

u 



undamped min 



1 



W cm -d~Ef> ( t 



/;, 



24(/i-3ff„) 

J dampedave 



v x - , 



(146) 



undamped max 



24( j u + 3<t |1 ) 



( t 


j 


2 


L 



(147) 



where <J fl is the estimated one sigma uncertainty of the 
value of fx. With 99.7% confidence, assuming a 
Gaussian distribution (and assuming u w(f) , u 



w(f>> u w u p(f)> 



■hi/p 



, and u E( f) equal zero), the value for f u 



U WC 1 J undampedmin tl *l" J imdampedmax ' 

calculated as follows: 



lies 
Therefore, u^, (f) is 



undamped 



/. 



undamped i 



J Hi 



(148) 



undamped i 



(the recommended approach), consult Sec. 3.2. This 
current section presents the equations for a step height 
measurement from two step height test structures on the 
same chip. 

To determine a step height measurement from two 
step height test structures on the same chip, the follow- 
ing three calculations are made: (1) a calculation of the 
reference platform height for both step height test 
structures, which assumes that the reference platforms 
are comprised of the same layers, (2) a calculation of 
the two platform heights involved in the step, and (3) a 
calculation of the step height. 

The first calculation is of the reference platform 
height for both of the step height test structures. To 
calculate the reference platform height, such as platlr 
from test structure "1," the following equations are 
used: 



(platlrWa + platlrWb + platlrWc) 
platlrW = ^- '-, (149) 



(platlrEa + platlrEb + platlrEc ) 
platl rE = ^- '- , (1 50) 



and 



*mn 



platlr - 



(platlrW + platlrE) 



(151) 



9. Appendix C. 

Supplemental Material for Step 
Height Measurements 

This appendix contains supplemental material for 
step height measurements. In Sec. 9.1, equations are 
presented for step height measurements taken from two 
step height test structures. In Sec. 9.2, the combined 
standard uncertainty equation is presented for this 
measurement. And, Sec. 9.3 presents the derivation 
of u Lstep , a component in the combined standard 
uncertainty calculation for a step height measurement 
taken from one step height test structure, as indicated in 
Sec. 4.3. 

9.1 Step Height Measurements From Two Step 
Height Test Structures 

Step height measurements can be taken from either 
one or two step height test structures. To obtain a step 
height measurement from one step height test structure 



where platNxD, platNxDt, and platNX in the List of 
Symbols provide details associated with the nomen- 
clature used. Similarly, to calculate the reference 
platform height, such as platlr from test structure "2," 
the following equations are used: 

(platlrWa + platlrWb + platlrWc) 
platlrW = ^- '- , (152) 



platlrE ■ 



(platlrEa + platlrEb + platlrEc ) 



and 



platlr ■■ 



(platlrW + platlrE) 



(153) 



(154) 



The second calculation is of the two platform heights 
(namely, platNX and platMY) involved in the step. The 
following equations are used: 

(platNXa + platNXb + platNXc) 
platNX = ^- '- - platNr 

(155) 
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and 

(platMYa. + platMYb + platMYc) 
platMY = ^- '- - platMr , 

(156) 

where platNXt in the List of Symbols provides details 
associated with the nomenclature used. Therefore, for 
platform "A" in test structure "1" and platform "B" in 
test structure "2," the corresponding platform heights 
(namely, platlA and platTS) are calculated using the 
following equations: 

(platl Aa + platl Ab + platlAc ) 
platlA = platlx 



1 platNX 



(157) 



and 



(platlBa + platlBb + platlBc ) 

platlB = - platlx . 

3 

(158) 

The third calculation is of the step height, stepN x M Y , 
where the List of Symbols provides details associated 
with the nomenclature. The following equation is used: 



J 



11 LplatNX "*" U WplatNX + U cen ^ U repeat + U drtft "*" ' U linear ' 



(162) 

where u LplatNX , u WplatNX , u ca1 , w repeat , w^g, and Mi^ar a re 
given in the List of Symbols and discussed in the 
paragraphs that follow and where u plmMY in 
Eq. (1 61) is found by replacing occurrences of N with 
M and occurrences of X with Y, as appropriate, in 
Eq. (162) and other applicable calculations. In deter- 
mining the combined standard uncertainty, a Type B 
evaluation [17] (i.e., one that uses means other than the 
statistical Type A analysis) is used for each source of 
error, except where noted. 

The uncertainty equation for u LplatNX is determined 
from the minimum and maximum platform height 
values (namely, platNX miR and platNX m3X , respectively) 
as given below: 



and 



platNX mm = \platNX\- 2>a LplatNX (163) 



platNX mm = \platNX\ + 3o LplalNX , ^ 



stepN x M Y = platMY - platNX , 



(159) 



where LplatNX , as derived within Sec. 9.3, is given by: 



which for a step height from platform 'A" in test 
stucture "1" to platform "B" in test structure "2" (or 
stepl A 2 B ) becomes the following: 



stepl A 2 B = platlB — platlA . 



(160) 



" LplatNX VV plttNXme "rougWX / " platNXzvt; "roughiVX ' 

(165) 

where <J platNXme is calculated using the following 
equation: 



" platNXs ~T~®platNX\> ~*~ ® platNX c 



platNXave 



(166) 



The combined standard uncertainty equation for the 
above calculation is given in Sec. 9.2. 

9.2 Uncertainty Calculation When Two Step 
Height Test Structures Are Used in a 
Step Height Measurement 

If two step height test structures on the same chip 
are used to obtain a step height measurement, such 
as stepNyMy from platNX and platMY as given in 
Eq. (159) in Sec. 9.1, the following equation is used to 
calculate the combined standard uncertainty, w cSH [17]: 



:fi 



2 2 

platNX "*" U platMY 



(161) 



where u platNX is given in the List of Symbols and is 
found using the following equation: 



and where O mughNX is the smallest of all the calibrated 
values obtained for O platNXt . However, if the surfaces of 
platNX, platMY, platNr, and platMr all have identical 



compositions, then c, 



rougtLVA" 



equals cr, 



rougM/7 > 



which 



equals the smallest of all the values obtained for 

GpIatNXt* GplatMYt> GplatNrDtt an d GplatMtDf With 99.7 /o COnll- 

dence, assuming a Gaussian distribution (and assuming 

equal zero), the value 
for \platNX\ lies between platNX min and platNX mim . 
Therefore, u LlatNX \s calculated as follows: 



U WplatNx' M « rt ' M repeat> M drift> atl d Wjj 



P^tNX max -platNX w 



(167) 



^LplatNX 



which simplifies to the following equation: 



U LplatNX ^LplatNX ® platNXave ^roughNX • \\vo) 
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The uncertainty equation for u WplatNX is determined 
from the minimum and maximum platform height 
values (namely, platNX miR and platNX , respectively) 
as given below: 



platNX^ = \platNX\ - 3<7 



WplaINX 



(169) 



and 



platNX min = \platNX\-3\platNX\ 



platNX max = \platNX\ + 3\platNX\ 



cert 



cert 



(175) 



(176) 



and 



platNX mm = \platNX\ + 3o WplatNX , (170) 

where <J WplatNX \s given by the following equation: 

4 



' WplatNX 



2 2 

I * platNX "*" S platNr > 



(171) 



where s platNX is the standard deviation of the platform 
measurements (namely, platNXa, platNXo, and 
platNXc) from the data traces and s pUNt is the standard 
deviation of the measurements used in the calculation 
of the reference platform, platNr. With 99.7 % confi- 
dence, assuming a Gaussian distribution (and assuming 
ULpiatwc, "«■«, "repeat, "drift, and M linear equal zero), the value 
for \platNX\ lies between platNX min and platNX^. 
Therefore, u WplatNX \s calculated as follows: 



PlatNX m - platNX m 



'IVplalNX 



(172) 



which simplifies to the following equation: 



■{• 



U WplaINX ® WplatNX V S platNX ~*~ S platN r • (1 73) 



Therefore, the following equation is used to find 
UwpiaAK ioxplatlA: 



u 



WptatlA 



~®WplatlA \ S platlA ~^ S platlr ' (174) 



where s phaA is the standard deviation of the measure- 
mevAsplatl Aa,platl Ab, andplatl Ac and where s platlr is 
the standard deviation of the measurements (namely, 
platlrWa, platlrWb, platlrWc, platlrEa, platlrEb, and 
platlvEc) used to find platlx. This component, u WplatlA , 
can be considered a Type A component. 

The uncertainty equation for u cert , is determined from 
the minimum and maximum platform height values 
(namely, platNX miR and platNX mia , respectively) as 
given: 



The uncertainty of the platform height is assumed to 
scale linearly with height. With 99.7 % confidence, 
assuming a Gaussian distribution (and assuming u LplatNX , 
Uw P MNx, "repeat, "drift, and w llnear equal zero), the value for 
\platNX\ lies between platNX miR and platNX m!a . 
Therefore, u cert is calculated as follows: 



PlatNX max - platNX a 



which simplifies to the following equation: 

a,. 



cert 



'-\platNX\ 



(177) 



(178) 



The above equation is comparable to Eq. (39) after 
replacing \stepN XY \ with \platNX\. 

The uncertainty equation for w rep( , at is determined 
from the minimum and maximum platform height 
values (namely, platNX min and platNX mia , respectively) 
as given below: 



and 



platNX min = \platNX\ - \platNX\ 



platNX^ = \platNX\ + \platNX\ 



2z, 



2z, 



(179) 



(180) 



The uncertainty of the platform height is assumed to 
scale linearly with height. Assuming a uniform distri- 
bution (and assuming u 



u 



LplatNX) u WplatNX> u cert> "drift, 



It,, 



U,<, 



and 



"linear equal zero), the value for \platNX\ lies between 
platNX min and platNX mia . Therefore, w repeat is calculated 
as follows: 



platNX mm - platNX n 



(181) 



which simplifies to the following equation: 

U rep e a,=^j==r\platNX\. (lg2) 
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The above equation is comparable to Eq. (43) after 
replacing \stepNxy\ with \platNX\, 

The uncertainty equation for H drift is determined from 
the minimum and maximum platform height 
values (namely, platNX miR and platNX m!ix , respectively) 
as given below: 



which simplifies to the following equation: 



and 



platNX n 



platNX n 



I platNX I - \platNX I -^f 

2z 



\platNX\ + \platNX\ 



Tz 



(183) 



(184) 



The uncertainty of the platform height is assumed to 
scale linearly with height. Assuming a uniform distri- 
bution (and assuming u LplatNX , u WpkaNX , u cert , u apedt , and 
"linear equal zero), the value for \platNX\ lies between 



platNX mm and platNX m 
follows: 



. Therefore, w drift is calculated as 



_ platNX max - platNX min 
" dnft= ~^3 ^' 

which simplifies to the following equation: 



\platNX \. 



lS~2 



(185) 



(186) 



The above equation is comparable to Eq. (47) after 
replacing \stepNxy\ with \platNX\. 

The uncertainty equation for w linear is determined from 
the minimum and maximum platform height values 
(namely, platNX miR and platNX mia , respectively) as 
given below: 



V^ 



\platNX\ 



(190) 



The above equation is comparable to Eq. (50) after 
replacing \stepNxr\ with \platNX\. 

9.3 Derivation of u Lstep 

As mentioned in Sec. 4.3, one component of the 
combined standard uncertainty equation for a step 
height measurement taken from one step height test 
structure is u Lstep , which is the uncertainty of the meas- 
urement across the length of the step. The derivation of 



■ l Lstep 



in Eq. (35) is given in this section. 



Recall that the step height as obtained from one step 
height test structure is given by Eq. (12) and repeated 
below: 



stepN x 



stepN^ + stepN xrb + stepN x 



(191) 



This equation can be written in terms of the platform 
measurements by using Eq. (8) to obtain the following 
equation: 

stepNx = -\(platNYa - platNXa) 

3 LV ^ F (192) 

+{platNYb - platNXb) +(platNYc -platNXc)] . 



Rearranging, this equation becomes: 

stepNxy = -[(platNYa. + platNYb + platNYc) 

3 (193) 

-(platNXa + platNXb + platNXc)] . 



platNX mm = \platNX\-\platNX\z pi 



and 



platNX^ 



\platNX\ + \platNX\ z 



(187) 



(188) 



The uncertainty of the platform height is assumed to 
scale linearly with height. Assuming a uniform distri- 
bution (and assuming u 



u 



LplatNXi u WpIatNX> u certi "repeat! 



U,., 



II,, 



and 

w drift equal zero), the value for \platNX\ lies between 
platNX min and platNX max . Therefore, w linear is calculated 
as follows: 



platNX ma - platNX m 



(189) 



Defining platNX me and platNY me as the average of the 
platform height measurements from platNX and 
platNY , respectively, this equation can be rewritten as 
follows: 



stepN ' M = platNY mc - platNX^ 



(194) 



Given the above equation, the uncertainty equation for 
ULstep can be written as follows: 



U L s,ep = V" 



LpNX ~*~ U LpNY > 



(195) 



where u LpNX is the uncertainty across the length of 
platNX and where u LpNY is the uncertainty across the 
length of platNY. 
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The uncertainty equation for u LpNX is determined from 
estimated minimum and maximum step height values 
(namely, stepNxy m{n and stepN^^, respectively) as 
given below: 



and 



stepN x 



stepN x 



platNY wc -[platNX. ive 



" ■* V platNXave "roughM" )i 

-- platNY me ~[platNX me 



-3((7 



platNXave 



-o„ 



ighM 



)• 



(196) 



(197) 



where o platNXmt and a r0llghA , z are defined in Sec. 4.3 and 
also in Sec. 9.2. With 99.7 % confidence, assuming a 
Gaussian distribution (and assuming u LpNY and all other 
sources of uncertainty equal zero), the value for stepNxy 
lies between stepNxy min and stepNxy^ . Therefore, u LpNX 
is calculated as follows: 



Ste P N XYm^- Ste P N X 



l LpNX 



(198) 



which [using Eq. (196) and Eq. (197)] simplifies to 



U LpNX ® platNXsve "roughWX ' 



(199) 



which is comparable to Eq. (165) in Sec. 9.2 for one 
platform. Similarly, the equation for u LpNY is as 
follows: 



U LpNY ^plaNYmc "roughAT 



(200) 



Inserting Eq. (199) and Eq. (200) into Eq. (195) 
produces the following equation for u Lstep : 



U Lslep V [yplalNXmc "roughWJf / "*" \y platNYswc "roughAT / 



(201) 



which is Eq. (35) in Sec. 4.3. This model for u Utep was 
chosen because if a plalNXme = a rm0MX = and a pIatNYme = 
"mughATi then u Lstep = 0. 
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